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Abstract 

Zeros of Random Analytic Functions 
by 

Manjunath Krishnapur 
Doctor of Philosophy in Statistics 

University of California, Berkeley 

Professor Yuval Peres, Chair 

The dominant theme of this thesis is that random matrix valued analytic functions, 
generalizing both random matrices and random analytic functions, for many pur- 
poses can (and perhaps should) be effectively studied in that level of generality. 

We study zeros of random analytic functions in one complex variable. It 
is known that there is a one parameter family of Gaussian analytic functions with 
zero sets that are stationary in each of the three symmetric spaces, namely the 
plane, the sphere and the unit disk, under the corresponding group of isometries. 

We show a way to generate non Gaussian random analytic functions 
whose zero sets are also stationary in the same domains. There are particular 
cases where the exact distribution of the zero set turns out to belong to an im- 
portant class of point processes known as determinantal point processes. 

Apart from questions regarding the exact distribution of zero sets, we also 
study certain asymptotic properties. We show asymptotic normality for smooth 
statistics applied to zeros of these random analytic functions. Lastly, we present 
some results on certain large deviation problems for the zeros of the planar and 
hyperbolic Gaussian analytic functions. 



Professor Yuval Peres 
Dissertation Committee Chair 
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Chapter 1 




1.1 Introductory remarks 

Random analytic functions on the one hand and random matrices on the 
other are two well studied topics in probability theory and mathematical physics. 
One of the chief interests to a probabilist in these objects is the kind of point 
processes one gets, by taking the set of zeros or the set of eigenvalues, as the 
case may be. Both these kinds of point processes typically have the property 
of "repulsion", meaning that the points distribute themselves more evenly than 
they would if they were thrown down independently. That is an appealing feature 
because, while there are ways to construct point processes that are more clumped 
than independent points, there are not many natural ways in which to generate 
point processes with less clumping. 

This fact and several others (more empirical than mathematical) have led 
to a folk wisdom that random analytic functions and random matrices share many 
similarities. Differing responses to this statement have been heard, including one 
that points out the obvious tautology here (after all the characteristic polynomial 
of a random matrix is a random polynomial), and another that says that there is 
not much similarity but instead evokes an "anthropic" reasoning (the same set 
of people work on both these fields). Without denying the validity of these expla- 
nations, in this thesis we take a more positive approach attempting to provide 
a unifying framework that includes both random matrices and random analytic 
functions (see caveats below). 
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This is the simple but seemingly useful idea of considering random matrix-valued 
analytic functions, and the set of points where it becomes singular (i.e., the zeros of 
the determinant). The linear polynomials reduce to random matrices and the lxl 
matrices correspond to random analytic functions. 

A little explanation is in order. When we talk of random analytic func- 
tions, we tacitly mean that we are somehow specifying the distribution of coef- 
ficients or some closely related quantities (otherwise any random set of points 
would be the zeros of a random analytic function) . Furthermore it is usually dif- 
ficult to analyse a random analytic function (especially to get exact properties) 
except in the case of Gaussian coefficients. So the essence of the above para- 
graph is that the determinant of a Gaussian matrix-valued analytic function is a 
non-Gaussian analytic function in itself, but nevertheless amenable to analysis 
because it is built out of Gaussian analytic functions. 

Secondly, the earlier claim about random matrices falling within our 
framework should be toned down. The chief, although not the whole, empha- 
sis in random matrix theory is on the study of Hermitian random matrices and 
their (real) eigenvalues, for physical as well as mathematical reasons. When we 
go to higher polynomials there is perhaps no natural way to get the zeros to lie on 
the real line. This may explain why these objects have not been studied before. 
What we study here are zeros in the complex plane, for which of course there is no 
such problem. Nevertheless we believe that it is also interesting mathematically 
to study polynomials with random Hermitian or random unitary coefficients (we 
do not do this here) even though the zeros are spread out in the complex plane. 

We now outline the contents of the thesis briefly. 

• In the remaining sections of this chapter, we give a quick introduction to 
the basic notions of a point process, correlation functions and Gaussian 
analytic functions. Most importantly, we recall the three canonical families 
of Gaussian analytic functions on the plane, the sphere and the unit disk 
(hyperbolic plane). 

• In Chapter|2]we give a recipe for generating a slew of (non-Gaussian) random 
analytic functions whose zeros are stationary in the plane, the sphere and 
the unit disk. We make some basic computations on the distribution of zeros 



3 



that will be used later. 

• In Chapter |3] we recall the notion of a determinantal point process, and 
characterize the stationary determinantal point processes in the three fun- 
damental domains. Of these the planar ones are known to be (limits of) the 
distribution of eigenvalues of certain random matrices (the Ginibre ensem- 
ble) while the processes on the sphere and disk are new (these processes 
themselves have been considered before in but an independent proba- 
bilistic meaning was not known) . 

• In Chapter |4[ we present an evocative analogy which suggests that the de- 
terminantal point processes on the sphere and the disk, introduced in Chap- 
ter^ are in fact the singular points of certain random matrix-valued analytic 
functions that were introduced in Chapter |2J 

• In Chapter [5] we prove that the stationary determinantal processes on the 
sphere introduced in Chapter |3] are the singular points of the random matrix 
analytic function zA — B (in this case, more simply, the eigenvalues of A~ 1 B). 

• In Chapter [6] we give partial proof that the determinantal processes on the 
disk introduced in Chapter |3] are the singular points of the random matrix 
analytic function A + zA\ + z 2 A 2 + 

• In Chapter [7] we show asymptotic normality for smooth statistics applied to 
the zeros of random analytic functions introduced in Chapter |2] following 
a method of Sodin and Tsirelson who showed the same for the canonical 
models of Gaussian analytic functions. 

• In Chapter |8] and Chapter |9] we move away from the line of presentation 
so far, and return to canonical Gaussian analytic functions. We deal with 
two large deviation type problems for zeros of the planar Gaussian analytic 
function, one posed by Yuval Peres, which we solve fully and another due to 
Mikhail Sodin, which we solve partially. 
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1.2 Basic notions and definitions 

1.2.1 Point processes, Correlation functions 

A point process in a locally compact Polish space 17 is a random integer- 
valued positive Radon measure X on Vt. (Recall that a Radon measure is a Borel 
measure which is finite on compact sets.) If X almost surely assigns at most mea- 
sure 1 to singletons, it is a simple point process; in this case X can be identified 
with a random discrete subset of f2, and X(D) represents the number of points of 
this set that fall in D. 

The distribution of a point process can, in most cases, be described by 
its correlation functions (also known as joint intensities) w.r.t a fixed Radon 
measure /i on Q. 

Definition 1.2.1. The correlation functions of a point process X w.r.t. fi are func- 
tions (if any exist) p n : Vt n — >■ [0, oo) for n > 1, such that for any family of mutually 
disjoint Borel subsets D\, . . . , of O, and for any non-negative integers n%, . . . , nk, 



where 



n 



E 



k 

i=l 



n 

i=l 



X(Di 



p n (x\, . . . , x n )dfj,(xi) . . . dn(x n ), 



(1.2.1) 



n 



Remark 1.2.2. It is a natural question to ask for conditions that guarantee the 
existence of correlation functions and conditions under which they determine 
the dist ribution of the point process. Such conditions do exist, see Lenard's 



23),I2 



25l) or the survey by Soshnikov 



3711 . But the conditions are too com- 



plicated and not relevant for our purposes. In any case, when the joint distribu- 
tion of X(Di), . . . ,X(Dk) is determined by its moments, the correlation functions 
determine the distribution of X. 



Remark 1.2.3. For overlapping sets, the situation is more complicated. Restrict- 
ing attention to simple point processes, p n is not the intensity measure of X n , but 
that of X An , the set of ordered n-tuples of distinct points of X. Indeed, 11.2. II 
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implies (see |23j; |24; 28)) that for any Borel set BcO"we have 

E#(BHX An ) = J p n (x u ...,x n )dn(xi)...d f i(x n ). (1.2.2) 

B 

Assuming that X is simple, the correlation functions may be interpreted 
as follows: 

• If 17 is finite and \i = counting measure then pk(x\, . . . ,x^) is the probability 
that xi, . . . , Xk € X. 

• If fi is open in M. d and p = Lebesgue measure, if p n exist and are continuous, 
then 

P\X has a point in each of BJxj)] , 

»(«..-.»*)- is- 1 i^p# — ■ (1 - 2 - 3 ' 

Conversely, if for every k > 1, the right hand side of 11.2.31 exists and is 
continuous in X{, 1 < i < k, then it is the /c-point correlation functions of X. 

For us f2 will always be an open subset of the plane (or the sphere S 2 ) 
and X will be a simple point process, pt may always be taken to be the Lebesgue 
measure on f2, but we often find it convenient to use some other measure that is 
mutually absolutely continuous with the Lebesgue measure. 



1.2.2 Complex Gaussian distribution 

A standard complex Gaussian is a complex- valued random variable with 
probability density ^e~' 2 ' 2 w.r.t the Lebesgue measure on the complex plane. 
Equivalently, one may define it as X + iY, where X and Y are i.i.d. N(0, \) random 
variables. 

Let a,k, 1 < k < n be i.i.d. standard complex Gaussians. Let a denote 
the column vector (ai, . . . , a n )*. Then if B is an m x n matrix, Ba + \i is said to be 
an m- dimensional complex Gaussian vector with mean fi (an m x 1 vector) and 
covariance E = BB* (an m x m matrix). We denote its distribution by CN m (fj,, S). 

Here are some basic properties of complex Gaussian random variables. 

• If a is a complex Gaussian, its distribution is determined by \x = E [a] and 

E = E [(a — fj,) (a — //)*]. All moments of the form 

E[(a fe - Hk){aj - Hj)] , l<k,j<n, 
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vanish. This is the case even for j = k. 

If a is a standard complex Gaussian, then \a\ 2 and A are independent, and 
have exponential distribution with mean 1 and uniform distribution on the 
circle {z : \z\ = 1}, respectively. 

If a n , n > 1 are i.i.d. CN(0, 1), then by an easy application of Borel-Cantelli, 
lim sup \a n \" = 1, almost surely. (1-2.4) 

n— >oo 

In fact, equation 11.2.41 is valid for any i.i.d. sequence of complex-valued 
random variables a n , such that E [max {log |oi|, 0}] < oo. Equation 11.2.41 
is useful to compute the radii of convergence of random power series with 
independent coefficients. 

Wick Expansions: The Wick or the Feynman diagram expansion is an ex- 
pansion of I? functions of a Gaussian measure in an orthonormal basis 
consisting of polynomials of the u nde rlying Gaussians. Following the pre- 
sentation in the book by Janson , we state the essential facts in the 
limited context that we shall need later. More details and complete proofs of 
the assertions can be found in 

Let ax,...,ap be i.i.d. standard complex Gaussians. Consider the collection 
of all monomials FT a™ k a% k in these p variables, and orthonormalise them by 

k 

projecting the polynomials of degree N on the orthogonal complement of the 
polynomials of degree N — 1. 

This procedure is the same as applying Gram-Schmidt to the monomials 
after arranging them in increasing order of the degree (how we order mono- 
mials of the same degree is immaterial because distinct monomials of the 
same degree are clearly orthogonal) . Thus we get an orthonormal basis of all 
square integrable functions of the a^s, and the basis elements, termed Wick 
powers, are denoted by 

mi. — ni. mt-tii. 
a u CLu T-r a k a k 



■ TT ? fe =TT 

11 11 



,, y/m k \n k \ L £ ^/ mfc !n fe ! 
the equality a consequence of the independence of a^s. 
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These Wick polynomials are known explicitly (see 117))- 



ml\n 





(-l)V! 



although this is not particularly important to us. It is quite well known that 
products of random variables that are jointly Gaussian can be described 
by summing over the weights of certain combinatorial entities. There is a 
similar formula (known as Wick formula or Feynman diagram formula) for 
expectation of product of Wick powers. We shall only need the following 
special case. 

Wick/Feynman diagram formula: Let (6j, . . . , b s ) have a complex Gaussian 
distribution with mean zero. Then 



where the sum is over all complete Feynman diagrams 7 without self inter- 
action (henceforth we shall just say Feynman diagram). To define this, 
consider a collection of XX m j + n j) vertices with rrij of the vertices labeled 



j and nj of the vertices labeled j, for 1 < j < s. All the vertices labeled j 
are also supposed to be distinguishable although we shall not introduce any 
more notation to distinguish them. Now, each 7 is a matching of these ver- 
tices (a subgraph in which each vertex has degree 1), such that each edge in 
7 connects a vertex labeled i to a vertex labeled j for some i / j. 

The value v(-f) of the diagram is the product of the weights of all the edges 
in 7, and the weight of an edge joining a vertex labeled i to a vertex labeled j 
[i + j) is E [bibj] . 

Example 1.2.4. Let s = 2. Then we must consider Feynman diagrams on 
the labels {1, 1,2,2} with mi vertices labeled 1, m,2 vertices labeled 2, ri\ ver- 
tices labeled 1 and 712 vertices labeled 2. Since a Feynman diagram (in our 
terminology as explained above) must connect Is to 2s and vice-versa, and 
must give every vertex degree one, there are no Feynman diagrams unless 



s 



E 




(1.2.5) 
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mi = n 2 and m 2 = n\, in which case there are m\\n\\ such diagrams. Thus 



E 



n : <w 



mi!ni!E[6i6 2 ] mi E[6 2 &i] ni if mi = n 2 ,m 2 = n\. 
otherwise. 



1.2.3 Gaussian analytic functions 

Endow the space of analytic functions on a region Q with the topology of 
uniform convergence on compact sets. This makes it a complete separable metric 
space which is the standard setting for doing probability theory (To see complete- 
ness, if {/„} is a Cauchy sequence, then /„ converges uniformly on compact sets 
to some continuous function /. Then it is easy to see that / must be analytic 
because its integral on any closed contour is zero since f f = lim f f n and the 

latter vanishes for every n, by analyticity of f n ). 

Definition 1.2.5. Let f be a random variable taking values in the space of analytic 
functions on a region c C. We say f is a Gaussian analytic function (GAF) on 
if (f(zi),... ,f(z n )) has a mean zero complex Gaussian distribution for every 

Z\ , . . . , z n G Q . 

It is easy to see the following properties of GAFs. 

• {fW} are jointly Gaussian, i.e., the joint distribution of f and finitely many 
derivatives of f at finitely many points, 

{f (fe) (^) :0</c<ra, l<j< m j , 
has a (mean zero) complex Gaussian distribution. 

• The distribution of a Gaussian analytic function is determined by its covari- 



ance kernel ( E 



i{z)i{w) 



denoted by K t (z,w) or just K(z,w) if there is 



no ambiguity as to which f is under consideration. 



1.2.4 Stationary zero sets of Gaussian analytic functions 

Our interest is in the zero set of a random analytic function. Unless 
one's intention is to model a particular physical phenomenon by a point process, 
there is one criterion that makes some point processes more interesting than 
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others, namely, stationarity under a large group of transformations (stationarity 
of a random process means invariance of its distribution under a group action. 
It is also called invariance, especially when the stationarity is in "space" rather 
than "time", but we use both terms interchangeably). There are three particular 
two dimensional domains on which the group of conformal automorphisms act 
transitively (There are two others that we do not consider here, the cylinder or the 
punctured plane, and the two dimensional torus). We introduce these domains 
now. 

• The Complex Plane C: The group of transformations 

<p x ^{z) = \z + (3, zgC (1.2.6) 

where |A| = 1 and (3 G C, is nothing but the Euclidean motion group. These 
transformations preserve the Euclidean metric ds 2 = dx 2 +dy 2 and the Lebesgue 
measure dm(z) = dxdy on the plane. 

• The Sphere § 2 = C U {oo}: The group of transformations 

<PaA z )= _ > zG<CU{oc} (1.2.7) 

— pz + a 

where a, f3 G C and |a| 2 + 1/?| 2 = 1, is the group of linear fractional transforma- 
tions mapping C U {oo} to itself bijectively. These transformations preserve 
the spherical metric ds 2 = p^n^pfp and the spherical area measure n^rfpL ■ 
We call it the spherical metric because it is the push forward of the usual 
metric on the sphere inherited from M. 3 , onto Cu{oo} under the stereographic 
projection, and the measure is the push forward of the spherical area mea- 
sure. The transformations 11.2.71 are just the rotations of the sphere under 
this identification with C U {oo}. 

• The Hyperbolic Plane B = {z : \z\ < 1}: The group of transformations 

i p a ^(z) = ^±^-, zeB (1.2.8) 

pz + a 

where a, (3 G C and \a\ 2 — \(3\ 2 = 1, is the group of linear fractional transfor- 
mations mapping the unit disk B = {,z:|2:|<l}to itself bijectively. These 
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transformations preserve the hyperbolic metric ds 2 = ^rcpp and the hyper- 
bolic area measure nrrr^pp (this normalization differs from the usual one, 
with curvature —1, by a factor of 4, but it makes the analogy with the other 
two cases more formally similar) . This is one of the many models discovered 
by Poincare for the hyperbolic geometry of Bolyai, Gauss and Lobachevsky 
(see 16) for an introduction). 

Note that in each case, the group of transformations acts transitively on 
the corresponding space, i.e., for every z,w in the domain, there is a transforma- 
tion ip such that (f(z) = w. This means that in these spaces every point is just 
like every other point. Now we introduce three families of GAFs whose relation to 
these symmetric spaces will be made clear in Proposition I1.2.7L 

In each case, the domain of the random analytic function can be found 
from 11.2.41 . Indeed, 11.2.41 implies that when a n are i.i.d. standard complex 
Gaussians, ^ a n c n z n has the same radius of convergence as ^ c n z n . 



• The Complex Plane C: Define for L > 0, 

:'L" 

n=0 

For every L > 0, this is a random analytic function in the entire plane. 

• The Sphere § 2 : Define for L e N = {1, 2, 3, . . .}, 



\/T n 

f(z) = Va n ^/. (1.2.9) 



f (*) = £ q^- 1 ) + (1.2. 10) 

— ' - '"7.! 



n=0 v'n! 



For every LeN, this is a random analytic function on S 2 = C U oo with a pole 
at oo (i.e., it is a polynomial). 

• The Hyperbolic Plane D: Define for L > 0, 



^^ ^MU.) ^ (1211) 

n=0 



For every L > 0, this is a random analytic function in the unit disk ID = {z : 

\z\ < 1}. 
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Remark 1.2.6. Although we wrote 11.2.91 for every L > 0, they are identical up to 
a scaling of the complex plane. However, the functions in 11.2. 101 and 11.2. 1 II are 
truly different for different L, i.e., there is no transformation of the S 2 and D, that 
makes and f^/ the same, for L / L'. This is particularly obvious for the sphere, 
because L then denotes the number of zeros of f . 



We just quote the following proposition from 13511 . (The proof is contained 
in the proof of Proposition 12. 1. II . These random analytic functions were discov- 
ered in several stages and (partially) by several authors. The main contributions 
are due to Bogomolny, Bohigas and Leboeuf II) and Q), Kostlan I22I . Shub and 



Smale 



3211 . Some of them are natural generalizations (to complex coefficients) of 



ran dom polynomials studied by Mark Kac in his founding papers starting with 



19). The special case L = 2, in the unit disk was derived also by Diaconis and 
Evans Q) as the limit of the logarithmic derivative of characteristic polynomials of 
random unitary matrices. The uniqueness in Proposition 1 1 . 2 . 71 also was perhaps 
known, but a much stronger form of uniqueness (that the first intensity of zeros of 
any Gaussian analytic function determines the distribution of the Gaussian ana- 
lytic function itself, up to multiplication by arbitrary deterministic non-vanishing 
analytic functions) was found by Sodin 33). 



Proposition 1.2.7. The zero sets of the GAF f in equations 11.2.91 . 11.2.101 and 
11.2.1 It are invariant (in distribution) under the transformations defined in equa- 
tions 1 1.2.61 , 11.2.71 and 11.2.81 respectively. This holds for every allowed value of 
the parameter L, namely L > 0/or the plane and the disk and L e N/or the sphere. 

Moreover, these are the only Gaussian analytic functions (up to multiplica- 
tion by deterministic non vanishing analytic functions) with stationary zero sets in 
these domains. 
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Chapter 2 

Stationary zero sets of random 
analytic functions 

As we saw in Proposition 1 1 .2.71 on each of the three domains C/S 2 /B, 
there is a one parameter family of Gaussian analytic functions whose zero sets 
are stationary under the corresponding group of isometries. Moreover, these are 
the only Gaussi an analytic functions on these domains with these properties. 
Indeed Hannay likens the uniqueness of the Gaussian analytic function in 
11.2.91 to that of the Poisson process or the thermal blackbody radiation. 

Here we stick to the three domains C/S 2 /B and ask for random analytic 
functions whose zero sets are stationary. By Proposition I1.2.7L we must nec- 
essarily seek among non-Gaussian analytic functions. A natural idea might by 
to replace i.i.d. Gaussians in the coefficients by i.i.d. complex-valued random 
variables from some other distribution. However, these seem difficult to anal- 
yse. Gaussian analytic functions have the nice property that the evaluations of 
the function and its derivatives are all Gaussian with distributions that we can 
explicitly work with and this fails in other cases. In fact we do not know of 
another example of a power series with i.i.d. coefficients whose zero set is sta- 
tionary (on any of these three domains) . We resolve this deadlock by constructing 
non-Gaussian analytic functions using Gaussian analytic functions as building 
blocks. 
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2.1 A recipe for stationary zero sets of random analytic 
functions 

Let Q be a (non-random) homogeneous polynomial in k variables with 
complex coefficients and let f be any Gaussian analytic function (not necessarily 
one of the canonical models denned in Section 11.2.41 . Then if fj,i < k are i.i.d. 
copies of f, then Q(fi, . . . ,fj.) is a random analytic function on the same domain 
as f . 

Proposition 2.1.1. Let Q be a homogeneous polynomial of degree d ink variables 
with complex coefficients, and let f be one of the canonical models of Gaussian 
functions in U.2.91 . \1.2.1Qt or M.2.1U . lfii,\<i<kare i.i.d. copies off, then the 
zero set of the random analytic function 

F(z) :=Q &(*),... ,f fc (2)) 
is stationary under the same group of isometries as the zero setoff. 

Proof. First we recall the proof of invariance of the zero set of the Gaussian ana- 
lytic functions in 11.2.91 . 11.2. 101 and 11.2. 1 II .. Fix an isometry 92 of £1 (given in 
11.2.61 . 11.2.71 and 11.2.81 1. In each of the three cases, there is a deterministic 
non-vanishing function A^x such that 

f(z) = A VtL (z)f (<p(z)) , (2.1.1) 

where in fact 

' e Lz\n\LW domain = C. 
&<p,l{z) = \ v'{ z )~* domain = S 2 . 

ip'(z)~^ domain = D. 
Note that the equality in 12.1.11 is for the entire process, not just for a fixed z. 
Therefore, the zero set of f is invariant in distribution under the action of (p. 
(To prove equation 12. 1. II . just compute the covariance kernels of the Gaussian 
processes on the left and right hand sides). 
Coming back to F, we see that 

A(z) d FMz)) = Q(A(z)f 1 (i P (z)),...,A(z)f k (i P (z))) 
= g(f x (*),..., ffcO)) from(|2XB 
= F(*). 
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This implies that the zero set of F is invariant in distribution under the action of 

ip. □ 

This is a very simple observation, but note that while F is built in a simple 
manner out of copies of f , the zero set of F is by no means a simple transformation 
of the zero sets of fi, . . . (except in trivial cases such as when Q(Ci, C2) = CiC2)- 
Thus the sets of zeros that we get are genuinely new point processes, but have the 
advantage of being based on Gaussian analytic functions, and therefore amenable 
to analysis. We illustrate this next, by computing the first and second correlations 
(joint intensities) for the zeros of F. The tool that we use to study functions such 
as F is the Wick expansion, suggested to us by Mikhail Sodin (see the paper by 



Sodin and Tsirelson 35) for a use of Wick expansions in the context of Gaussian 
analytic functions). We call random analytic functions of the kind described in 
Proposition 12. 1 . ll as polygafs. 



2.2 How to study the zeros of a polygaf? 

If F is any analytic function (not random) on f2, let dn F denote the count- 
ing measure, with appropriate multiplicities, on the zeros of F. Then, 

^-Alog|F(z)| =dn F {z) (2.2.1) 

in the sense of distributions. This just means that for any ip e C^°({1), 

J <p(z)dn F (z) = J A^(z)i log \F(z)\dm(z), (2.2.2) 
n a 

where m is the Lebesgue measure. Therefore when F is any random analytic 
function, understanding the distribution of the zero set depends on being able 
to do computations with log |F| (When F is Gaussian, there are other approaches 
to studying the zero set of F, but it appears that the approach outlined here is 
the only one that is equally convenient for our more general setting. The other 
methods make use of the probability density of F evaluated at several points in 
the domain etc, which are not available to us here). 
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Now from 12.2.21 . if cpi, 1 < i < k are smooth functions with disjoint sup- 
ports in Q, we get that 



E 



i=l 



nio g |F(^)i 

k 

nio g |F( 



i=i 



k 



i=l 



i=l 



In the last line we integrated by parts. 

In ll.2. ll we denned the correlation functions in terms of the moments of 
the joint counts of the number of points falling in several regions. Fixing k distinct 
points wi, . . . , Wk in ^ and letting tpi be a bump function in a small neighbourhood 
of Zi, by elementary measure theoretical arguments one can deduce that the k- 
point correlation function pk of the zero set of F, w.r.t Lebesgue measure is given 

by 



Pk(wi, ■ ■ ■ ,Wk) 



J]log|FK;)| 



.1=1 



(2.2.3) 



for distinct w\ , . . . , Wk . 

Remark 2.2.1. Hammersley gave a formula for the correlation functions of 
zeros of random polynomials in terms of the distribution of the coefficients. 12.2.31 
is an alternative way of expressing the same. In this form, for k = 1 it is sometimes 
called Edelman-Kostlan formula. 

The way to analyse log |F| is via Wick expansions that were outlined in 
Chapter □ 

Example 2.2.2. The particular example of Wick expansions that is of interest 
to us is the following: Let Q be a homogeneous polynomial in k variables with 
complex coefficients. If a,i,i < k are i.i.d. C/V(0, 1) random variables, then 



E 



log | Q(ai ,...,«*) 



< oo 



for every finite p. Hence we can expand log \ Q(ai, . . . ,afc)| in Wick powers as 



log|Q(ai,...,o fc ) 



\ ^ ( m .n "I r 

l.nPZ* .7 = 1 



: a j J aj 



(2.2.4) 
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where m = (mi, . . . , m^), n = (m, 



, nfc), m! = mj! and 



V m!n! 



E 



Iog|Q(oi,...,a fc )| JJ : a^a 

3=1 



mi „ n j . 
'3 ■ 



and the equality in 12.2.41 is in the L 2 sense (it could be better, of course) . 
We record two observations for later use. 

• C m>n = C nim for all m, n G 7L\, because a k are also i.i.d. CiV(0, 1). 

• Cm,n = unless m. = n., where m. := Yl m j- To see this, note for any A with 

3 

|A| = 1, Xaj are also i.i.d. CiV(0, 1) and hence, by the homogeneity of Q, it 
is also true that log \ Q(Xai, . . . , \cik)\ = log \ Q(ai, . . . , etfc)|. Therefore, from the 
equation above for C m n , we see that C m n 
true unless m, — n, or C m n — 0. 



i, which cannot be 



2.3 Distribution of the zero set of a polygaf 

Now let fj be i.i.d. copies of f , a Gaussian analytic function on a domain 
f2 (not necessarily one of the canonical GAFs on the plane, sphere or disk). As 
before Q is a homogeneous polynomial. 

Define F(z) = Q(fi(z), . . . , ffc(z)). If K(z,w) is the covariance kernel off, 

then set 

^ ) = K^|^ = Q ( fl(z) '---' fe 
where h(z) = . fj ^ — and d is the degree of Q. Then from 12.2.41 we can write, 



log|F(z)|= £ -^Ll{:f^)'"T ; (:)- : 
— , vm!n! 

m,ne^ J=l 

We work with F rather than F, because ij(z) are i.i.d standard Gaussians for any 
z, and so we can directly use the Wick formulas that we stated in the previous 
chapter. 

First Intensity: The first intensity (or 1 -point correlation function) as given by 
l2~2~3l is 

Pl (z) = -!-A 2 E[log|F(z)|]. (2.3.1) 

Z7T 
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Therefore, 



1 

1 

2^ 



A Z E [log|F(z 



1 



log |F(z)|J + — A,logiT(z,z) d / 2 



^A.Cco + ^A.log^z,^ 2 



from 12.2.41 . Therefore we obtain 

Pi (^) 



4"7T 



A 2 log 



(2.3.2) 



As a special case, set k = 1, Q(C) = C m formula 12.3.21 to deduce that the intensity 
of zeros of f is ^A z log K(z, z) (This is known as the Edelman-Kostlan formula). 
Thus we see that the intensity of zeros of F is d times the intensity of zeroes of f . 
This simple relationship between the intensities is surprisingly not quite obvious 
from the definition. 

Two point Correlations: Again from 12.2.31 we get the 2-point correlation. It is 
easy to see that 

1 



p 2 (z,w) - pi{z)pi(w) 



(27T) 



2 A 2 A„,E 



log|F(z)|log|F(» 



From 12.2.41 the right hand side can be written as 



(27T) 



2 A 2 A„, 



^ , Vmlnlm'ln'! 1J : L 3 3 3 3 

n,m',n' j=l 



m,n,m' ,n 



This is precisely the situation elucidated in Example II .2.41 Thus, only terms with 



Now make use of the observations made earlier- (1) C m n = C n>m , and (2) 
Cm,n = unless m. = n. . Grouping together terms by p = m. , we get, 

\K(z,w)\ 2p 



p 2 (z,w) - pi(z)pi(w) = j2~^^2\C p \ 2 A z A v 

^ 71 ' p=0 



K(z,z)pK{w,w)p' 



(2.3.3) 



where \C P \ 2 



|Cm,n| 

m.=n.=p 



Remark 2.3.1. Equation 12.3.31 has the appealing feature that the effects of the 
two ingredients of F, namely the polynomial Q and the Gaussian analytic function 



f, are clearly separated. \C P \ 2 depends only on Q while 



\K(z,w)\ 2 P 
K(z,z)PK(w,w)P 



depends only 
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on f . This observation is crucially used in the next chapter, when we compute 
the correlations for specific polynomials Q. One can write analogous but more 
complicated expressions for higher correlations, with A Zk , 1 < k < m applied to a 
sum over Feynman diagrams. 



19 

Chapter 3 

Stationary determinantal point 
processes 

In this chapter we move away from random analytic functions and talk 
about a different class of point processes. In the next chapter, we return to 
zeros of random analytic functions and show that there are point processes in the 
intersection of the two classes. 

One of the main qualitative properties of zero sets of random analytic 
functions is that they have the property of "repulsion", also called "negative cor- 
relation", at short ranges. This terminology is a little misleading because corre- 
lations are never negative! The precise meaning of negative correlations is that 
p2(x,y) < pi(x)pi(y) for x, y that are sufficiently close. There is another class of 
point processes that has this repulsion property at all distances in a very strong 
sense. These point processes were introduced by Macchi 126) and are known as 
Determinantal (Fermionic) point processes. See Figure |3] for a visual compari- 
son of zeros and eigenvalues with a Poisson process. 

Definition 3.0.2. A point process X on is said to be a determinantal process 

with kernel K : 0? —* C, if it is simple (i.e., there are no coincident points almost 
surely) and its correlation functions w.r.t a measure p satisfy: 

p k (xt, . . . ,x k ) = det (K^xj))^^ , (3.0.1) 

for every k > 1 and x±, . . . , x^ GO. We shall always assume that K is the projection 
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Figure 3.1: Samples of a translation invariant determinantal process(left) and 
zeros of a Gaussian analytic function(right) . Determinantal processes exhibit re- 
pulsion at all distances, and the zeros repel at short distances only. The picture 
in the center shows a Poisson process (picture due to Balint Virag) 



kernel of a closed subspace of L , i.e., that K(x, y) = Y1 <Pj( z )<Pj(,y) f° r a sequence of 

i 

functions (there may be infinitely many of them) {cpj} orthonormal in L 2 (fi). The 
distribution is determined by giving the kernel IK or the Hilbert space on which it 
is a projection. 

Remark 3.0.3. This definition may look artificial at first sight, but the motivation 
comes from quantum mechanics, where the probability densities are given by the 
absolute square of a complex- valued function called the amplitude. If ipi,...,<p n 
(assumed to be orthogonal) are single particle wave functions of electrons, the 
most natural re-particle wave function is not f\ <fik{xk) ( as would have been the case 

k 

for "independence") because firstly it has no symmetry in x±, . . . , x n , and secondly 
it shows no properties such as repulsion. Hence the idea is to anti-symmetrize 
this, to get det ((fi(xj)) (This is for Fermions. For Bosons, one symmetrizes and 
gets the permanent of (<p%(xj)) t but that is another story). The probability density is 
given by the absolute square of this, which can be written as det (K(xi,Xj)) 1<i j <n , 

n 

with K(x,y) = < -Pj{ x )Wj{y)- To generalize this notion to infinite particle systems, 
it is necessary to formulate the definition in terms of correlation functions as 
done in Definition 13.0.21 For correlation functions to be non-negative in 13.0. II a 
natural assumption is to take K to be a Hermitian non-negative definite kernel. 
Then it turns out that < K < 1 is a necessary condition for there to exist a 
determinantal process with kernel K (see I37I or Il4l ). And then, any such 
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process can b e expres sed as a mixture of determinantal processes with projection 
kernels (see 301. 13 ll) . ll4i ). The rank of the projection is the number of points 
in the process. Observe that p2(x,y) — p\(x)p\{y) = —\K(x,y)\ 2 is negative. More 
general inequalities like this are a clear consequence of the determinantal form of 
the correlation functions. 

Remark 3.0.4. Another justification for this definition is the huge number of 
instances of determinantal processes that have arisen so far in random matrix 
theory and combinatorics, many of them predating the definition. To name a few, 



Non-intersecting random walks on the line (Karlin and McGregor 



20), 



18)) 



Eigenvalues of a random unitary matrix chosen according to the Haar mea- 
sure (Dyson 19)). (There are many more random matrix examples). 

Subset of edges of a finite graph present in a uniformly chosen spanning tree 
(Burton and Pemantle 14)). 

(An encoding of) Young diagrams sampled from the poissonised Plancherel 
measure of the symmetric group (Borodin, Okounkov and Olshanski B). 



The extensive survey of Soshnikov 37) gives many more examples and details. 



While determinantal processes on M or Z have been studied extensively 
because they arise in random matrix theory and combinatorics, in two dimensions 
they seem to have been largely untouched. 

To get a determinantal point process is trivial. Just take a reproducing 
kernel Hilbert space H of functions on f2, and let K be the reproducing kernel of 
H. Then there exists a determinantal point process with kernel K (trivial when 
dim(H)< oo, otherwise it can be constructed by taking limits of finite dimensional 
ones. See 37) or for details). However to deserve our attention, the process 
must have some attractive features in addition to its determinantal nature. We 
adhere to two guiding principles. 

• That the process be invariant in distribution (i.e., stationary) under a rich 
group of transformations. 

• That the process arise in a natural way probabilistically. 
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Although imprecise, the first principle suggests that we consider the same three 
domains C,8 2 ,D. And since we would ultimately want to relate these to zeros of 
random analytic functions, we consider Hilbert spaces of analytic functions on 
these spaces. We show that the Hilbert spaces that give rise to stationary de- 
terminantal processes are precisely the well known Bargmann-Fock spaces of 
analytic functions. But before going into this, we should point out that these 
determinantal processes were already studied by Caillol 1 5) under the name "One 
component plasma on the sphere" (the two-component version was studied by 



Forrester, Jancovici and Madore 



1010 . and Jancovici and Tellez 16) from the point 



of view of constructing Coulomb gases on these spaces. We arrived at these pro- 
cesses prompted by a question of Balint Virag as to what natural determinantal 
processes can be defined on the two dimensional sphere. As we remarked ear- 
lier, one can define a determinantal process (which can be regarded as Coulomb 
gas at a particular temperature j3 = 2) by choosing one's favourite Hilbert space 
with a kernel. What is new here is that we show that these are unique in a cer- 
tain sense, and most importantly, we show in the next chapter how to get these 
determinantal processes as zeros of random analytic functions. 

Theorem 3.0.5. Let £1 be one ofC, § 2 , B and let fi be an arbitrary radially symmetric 
Radon measure on f2. Let ipj. be complex analytic functions on £1 and belong to L 2 (/i). 
Then the determinantal process with kernel 

N 

K{z, w)=J2 Vk^kH (3-0.2) 

k=l 

is invariant in distribution under the corresponding group of isometries if and only 
if it is one of the following. 

• $7 = C, dfi(z) = e~ a \ z \ 2 dm{z) , N = oo, ipk(z) = J~^z k , where a G (0,oo), and the 
kernel is 

K^(z,w) = -e azW . (3.0.3) 

7T 

The Hilbert space is the space of analytic functions in L 2 (C, ^e~ a l 2 l 2 ). We call 
this process Det-C-a . 



n = S 2 , dfi(z) = ( Ul , N = a, tp h (z) = J% z k , where a e {1,2,3,...}, and 

(l + H 2 )^ 
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the kernel is 

Kl'(z,w) = -{l + zw) a - 1 . (3.0.4) 

The Hilbert space is the space of analytic functions in L 2 (§ 2 , n ^i + ^\i)a+i ) • We 
call this process Det-S 2 -a . 

• n = D, dfi(z) = (1 - \z\ 2 )^ a -^dm{z), N = oo, <p fc (z) = ^/fz*. where a G (0,oo), 
and the kernel is 

The Hilbert space is the space of analytic functions in L 2 (D, ^(1 — |z| 2 ) a_1 ). We 
call this process Det-B-a . 

Remark 3.0.6. Note the similarity to the classification of Gaussian analytic func- 
tion with stationary zeros in Proposition 1 1 . 2 . 71 Just as there, here too, Det-C — a 
processes are all identical up to scale, whereas the Det-§ 2 — a and Det-B — a are 
genuine one parameter families of point processes. 

There are (at least) two ways of using a positive definite kernel in prob- 
ability theory. One is to use it as the covariance kernel of a Gaussian process 
and another is to use it as the kernel of a determinantal process (it has to be a 
projection kernel for the latter). We are not aware of any probabilistic connection 
between the two. 

When we use the reproducing kernels of the Bargmann-Fock spaces as 
covariance kernels of Gaussian processes, we get the Gaussian analytic functions 
introduced in 11.2.91 . 11.2. 101 and 11.2. 1 II . When we use them as kernels for 
determinantal processes we get stationary point processes in these domains. 

Proof ofTheorem &ZTEl Let 

i if = C. 



— /-, . I — TTY7 tf — § 2 . 



i 



if n 



\ 7r(l-|2| 2 ) 2 

Then p 2 (z)dm(z) is the unique invariant measure (up to multiplication by a con- 
stant) on 0. 

If X is a determinantal process with kernel I on n, with distribution 
invariant under the corresponding isometry group, then the first intensity of the 
process, K(z, z)dfi(z) must be equal to ap 2 (z)dm{z) for some a > 0. 
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Express the correlation functions of X w.r.t the measure p 2 (z)dm(z) in- 
stead of the Lebesgue measure. Then the kernel becomes 

aK(z, w) 

y/K(z,z)y/K(w,w)' 

Invariance of the second correlation function implies that 

K(<p(z), tp(z))K(<p(w), tp(w)) 

for every isometry tp of 0, where by Const(z, w) we mean that it does not depend 
on tp. 

The idea is this. We differentiate equation 13.0.61 w.r.t tp and equate to 
zero. The derivatives w.r.t tp can be written as derivatives w.r.t z, w. That gives us 
differential equations for K that are easy to solve. 

Firstly, fix any (zq,wq) such that K(zq,wo) / 0. Without loss of generality 
take z = = wq. Then there is a neighbourhood N of in the Complex plane 
such that if z,w G N, and tp is close to identity, then \K(z, w) - K(0, 0)| < ||K(0, 0)|. 
Let S be the disk of radius ||K(0, 0)| centered at K(0, 0). Then S cannot intersect 
both the positive and negative parts of real axis because S is convex and does not 
contain 0. Moreover S and S intersect the real line at the same points. Therefore 
by removing the positive or the negative half line, we can define a continuous 
branch of logarithm on S U S. Henceforth "log" will denote this function. 

Taking logarithms in Equation l3.0.6l we get 

logK(p(z),p{w)) +logK(p(w),p(z)) -logK(p{z),p(z)) -logK(p(w),p(w)) (3.0.7) 

is equal to const(z, w), not depending on tp. To differentiate w.r.t tp we parameterize 
it with complex numbers as follows. 



Complex plane Write tp(z) = Xz + a, where a £ C, |A| = 1. 



• Sphere Write tp{z) = where a,0 G C, \a\ 2 + |/3| 2 = 1. 

• Disk Write tp(z) = 2?±g, where a,P G C, Id 2 - I/3I 2 = 1. 

Let us deal with the planar case first. 
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Complex plane Write tp a (z) = z + a, where a G C. Then for e is small 
enough, if \a\ < e and z, w G iV, then K(<^ Q (z), ip a (w)) G 5. Apply g|^= to equation 
13.0.71 and evaluate at a = 0. We get 



2 



log K(z, w) + — — log K(w, z) - — — log K(z, z) - - log K(w, w) 



r\ r\ O \ ~ J ' r\ r\ — O \ 7 J t~\ r\ — O \ " J O O 

ozow uwuz ozoz owow 

= Q(z,w) + Q(w,z) - Q(z,z) - Q(w,w). 

Where Q(z,w) = ^= logK(z,to). This is well defined for z,w G N and is ana- 
lytic in z,w. Appling ^= to Q(z,w) + Q(w,z) — Q(z,z) — Q(w,w) we deduce that 
g^=Q(z,w) = 0. By expanding Q locally as a power series in z,w, and observing 
the symmetry Q(z,w) = Q(w,z), we conclude that Q(z,w) = h(z) + h(w) for some h 
holomorphic on N. Then 

log K(z,w) = whi (z) + zh\ (w) + C, \/z,w G N, 

where h[ = h and C is a constant. 

Now again consider 13.0.71 and apply ^ to it. We get 

d d d d 

= — \ogK{z,w) + — logK(w,z) - ^- logK(z,z) - — logK(w,w) 
oz ow oz ow 

= Wh(z) + hi(w) + zh(w) + hi(z) — zh(z) — h%(z) — wh(w) — hi(w) 
= (z-w)(h{w) -h(z)). 

Since this has to hold for every z,w G N, we must have h = 0. That means that 
hi(z) = az + b for some a, b. Therefore, 



\og~K(z,w) = w(az + b) + z{aw + b) + C 



= (a + a)(z + — *!—)( w + — h —) + C' 
a + a a + a 

Making a change of variables we get K(z,w) = e azw . Then a has to be positive. 
This concludes the planar case. 

Sphere In this case, we again differentiate Equation l3 . . 61 w.r.t a, (3 and 
their conjugates. However (p aj p depends on the parameters and their conjugates 
as well, and that makes the equation longer. A simplification is obtained by noting 
the following: 

Let g(z, w) be analytic in z,w. Then with ip = (p a ,p, 

(— +W— )(-= + z— )g(<p(z)Mw)) = - ( = R l_ v A bdtgXipWMw)). 
dp da d/3 da (—(3z + a)(—/3w + a) 

(3.0.8) 
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Here 81,82 denote the derivatives w.r.t the first and second arguments. 

Apply + w-j^)(jL + z-§=) to Equation 13.0.71 and evaluate at a = 1,(3 = 0. 
Then 13.0.81 yields, 

= Q(z, w) + Q(w, z) - Q(z, z) - Q(w, w), 

where, this time Q(z, w) = (1 + z w) 2 q^= logK(z,w). But all the considerations that 
applied to Q in the Planar case also apply here and we get K(z,w) = (1 + zw) a 
(perhaps after a change of variables). Now a will have to be a positive integer 
(because integrating K(z,z) over the whole space should give AH, where iV is the 
total number o points in the process) . 

Disk Analogous to the spherical case, here we observe that for any g 
analytic in z and anti-analytic in w, 

Q Q Q Q (\ — zw) 2 

(w- W to ){ W~ X fc )9MZ) '^ iW)) = (^ + a)(/^ + a) (ai ^ )( ^ (z) '^ H) - 

(3.0.9) 

Applying - w^)(M= - z -§=) to l3~07l and using equation <3^9l gives 

us 

= Q(z, w) + Q(w, z) - Q(z, z) - Q(w, w), 

where, this time Q{z,w) = (1 — zw) 2 -^=\ogK{z,w). As before this leads us to 
K(z, w) = (1 — zw) a and a will have to be positive. 

This completes the proof of the theorem. □ 
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Chapter 4 

Random matrix-valued analytic 
functions 

In Chapter |^1 we saw that by choosing a Gaussian analytic function f 
with a stationary zero set and a (non-random) homogeneous polynomial Q, we 
could construct a random analytic function with stationary zeros. There are two 
complementary questions that arise naturally. 

1. Can one study these random analytic functions in this generality without 
having to appeal to special Q and f ? 

2. Are there particular examples of Q and f that are somehow special? 

The answer to both these questions is yes. Regarding the first question, 
we already saw in Chapter |2] that the correlation functions of the zero set can be 
computed in a general fashion. We shall use these computations in Chapter [7] to 
prove asymptotic normality for the zero sets in general. In the current chapter 
and the next two, we answer the second question and show that zeros of random 
analytic functions sometimes (but far from frequently, let alone always) turn out 
to be determinantal point processes. 
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4.1 Determinantal processes that are zeros of RAFs: Known 
results 

We remarked earlier that the focus in random matrix theory has been 
on Hermitian random matrices. In the preface to his book "Random matrices", 



Mehta 27) says "The theory of non-Hermitian random matrices, though not ap- 
plicable to any physical problems, is a fascinating subject and must be studied 
for its own sake. In this direction an impressive step [has been taken by] Gini- 
bre ..." Ginibre found the exact distribution of eigenvalues of three (two, strictly 
speaking) ensembles of non-Hermitian random matrices. We quote the one that 
is relevant to us. 



Theorem 4.1.1 (Ginibre(1965) 



11)). Let Abe annxn matrix with ltd. standard 



complex Gaussian entries. Then the eigenvalues of A have density 

1 - E Nl 2 n , l2 
p n (zi, . . . ,z n ) = — — — e [JJzi - Zj| . (4.1.1) 

7T n J] fe! i<3 

k=l 

Equivalently, one may say that the eigenvalues of A form a determinantal point 
process with kernel 

n—l r n/% 

K n (z,w)=J2 [ -^T, (4.1.2) 

k=0 

w.r.t the reference measure dfi(z) = ^e - ' 2 ' 2 . The corresponding Hilbert space H = 
span{l, z,..., z^ 1 } C L 2 (C, £^-dm(z)). 

Despite the enthusiastic response, as shown by Mehta's quote above, 
there do not seem to be any significant e xac t results beyond Ginibre's. The fol- 
lowing beautiful result of Peres and Virag 28) seems to be the next such. 

Theorem 4.1.2 (Peres and Virag(2003) ^28^ ). Letf be the random power series 
whose coefficients are Ltd. standard complex Gaussians (this is the case £1 = 
D, L = 1 in 11.2. Hi . Then the zeros of f form a determinantal point process on the 
unit disk O with the kernel (the Bergman kernel of the unit disk) 

1 

K(z,w) = — 

7T(1 — ZW) Z 



29 



w.r.t the reference measure dfj,(z) = ^dm(z) on B. The corresponding Hilbert space 
H = span{l, z,z 2 . . . ,} C L 2 (B, dm ^ ) is the space of all analytic functions in L 2 (B). 

Remark 4.1.3. Observing that Theorem 14. 1 .21 identifies the distribution of zeros 
of the Gaussian analytic function (with L = 1) defined in 11.2. 1 II as being Det-B — 1 
(recall the definition of Det-B — 1 from 13.0.51 1. one is tempted to guess that the 
Gaussian analytic functions defined in 11.2.91 11.2. 101 and 11.2. 1 II might have 
zeros distributed like the Det-C - L, Det-S 2 - L and Det-B - L (defined in 13.0.31 . 
13.0.41 and !3.0.31 . respectively). However these canonical Gaussian analytic func- 
tions do not have determinantal zero sets. Indeed, it was observed by Peres and 
Virag in their paper that these zero sets do not have negative correlations at large 
distances and hence, cannot be determinantal). Therefore this beautiful might- 
have-been story is completely false! Nevertheless, to quote Einstein, "The Lord 
is subtle, but not malicious". In the next section we shall see how a completely 
different but equally compelling picture might well be true. 



4.2 Determinantal processes that are zeros of RAFs: An 
analogy 

First let us list all the Gaussian analytic functions whose zero sets we 
know to be determinantal. This includes Theorem 14. 1.21 and two trivial cases (a 
one-point point process is always determinantal!). 

• z — a = : One zero, with standard complex Gaussian distribution on 
C. H = span{l} in L 2 (C, ^f-dm(z)). 

• za — b = : One zero, distributed uniformly on S 2 upon stereographic 
projection from the plane. H = spanjl} in L 2 (§ 2 = C U l oo}, n ^ l _^ z ^ dm(z)) 



• a + za\ + z 2 a2 + • • • = : Peres and Virag 28): Infinitely many zeroes in 
the disk. A determinantal point process with kernel 7r _1 (l — |z| 2 )~ 2 w.r.t Lebesgue 
measure on the unit disk. Equivalently, H = spanjl, z, z 2 , . . .} in L 2 (B, -dm(z)). 

Note that Ginibre's result (Theorem 14. 1. II can be seen as regarding the 
zeros of the random analytic function zl — A, which can be thought of as a matrix 
version of the first of the above examples. This suggests that we consider the 
matrix versions of the other two, i.e., we look at 
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• det(zA-B) = 0, where A, B are n x n independent matrices with i.i.d. standard 
complex Gaussian entries. 

• det (^4o + zA\ + z 2 A2 + . . .) =0, where A\~ are independent n x n matrices with 
each one having i.i.d. complex Gaussian entries. 

The analogy strongly suggests that the solutions to these equations should 
give us the determinantal point processes corresponding to the Bargmann-Fock 
spaces on the sphere and the unit disk (but only for integer values of the pa- 
rameter, since the size of the matrix, namely n, is discrete). Before going into 
the statements and proofs, we make some big-picture remarks and connect these 
objects to the random analytic functions studied in Chapter |2j 

Remark 4.2.1. Note that here we are looking at the set of z for which a random 
matrix-valued analytic function (zA — B or Aq + zA\ + z 2 A2 + . . .) becomes singu- 
lar. This concept is an obvious generalization of both random matrices (which 
correspond to the case when the analytic function is linear) and Gaussian ana- 
lytic functions (which correspond to the case when the matrices have size lxl). 
In spite of this natural appeal, the concept of a random matrix-valued analytic 
function does not seem to have been considered in the literature. 

One possible reason could be that the focus in random matrix theory has 
been almost entirely on eigenvalues in one dimension (real line or the circle) for 
physical reasons as well as the strong mathematical connections with orthogonal 
polynomials, representation theory etc. Moreover the eigenvalues have a physical 
meaning in quantum mechanics. Note the difficulty of forcing the zeros to lie on 
the real line, except by considering eigenvalues of a Hermitian matrix. Never- 
theless, the idea of matrixifying seems to be useful, not only as suggested above 
with Gaussian matrix coefficients, but also polynomials with coefficients that are 
Haar -distributed unitary matrices. 

4.3 Matrix-valued GAFs and polygafs 

Now we want to point out the connection with homogeneous polynomials 
applied to i.i.d. copies of Gaussian analytic functions (polygafs, that is). 
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Consider det(zA — B). This is the same as applying the homogeneous 
polynomial Q ="det" in n 2 variables, to n 2 i.i.d. copies of the Gaussian analytic 
function {(z) = az - b (which is the case Q = § 2 , L = 1 in 11.2. 101 1. 

Similarly det [A + zA\ + z 2 A 2 + . . .) is the homogeneous polynomial Q =" 
det" in n 2 variables, applied to n 2 i.i.d. copies of the Gaussian analytic function 
f(z) = a + a\z + a 2 z 2 + . . . (which is the case f2 = B, L = 1 in 11.2. 1 II ). 

In other words, we have already shown in Proposition 12 . 1 . Tl that the zero 
sets of these RAFs are stationary in $7. In the next two chapters we investigate 
the distributions in greater depth. We shall show that in the first case ($7 = § 2 ) we 
do get determinantal processes, whereas in the second, we show partial results 
in this direction. The precise statements of the conjectures are as follows: 

Conjecture 4.3.1. Let A,B be ltd n x n matrices with ltd. standard complex 
Gaussian entries. The zeros of det(zA — B)form a determinantal point process with 
kernel 

n (1 + zw) n ~ l 



K(z,w) 



_ , , . n+1 . . ~ . n+1 1 

T(l + \z\ 2 ) — (1 + \w\ 2 ) — 



w.r.t the Lebesgue measure on C. Equivalently, K is the projection kernel on the 
subspace of analytic functions in L 2 (S 2 , n ^ 1+ ^^ n +i dm{zj). 

Conjecture 4.3.2. Let A/, be Ltd. nxn matrices with Ltd. standard complex Gaus- 
sian entries. The zeros of det (^4o + zAi + z 2 A2 + . . .) form a determinantal point 
process on B with kernel 

n (l - \z\ 2 ) 2 1 - \ w \*) i 
K{z,w) = — — , 

7T (1 — ZW) n+i 

w.r.t the Lebesgue measure on B. Equivalently, K is the projection kernel on the 
subspace of analytic functions in L 2 (B, - (1 — \z\ 2 ) n ~ 1 dm(z)). 

Remark 4.3.3. The Det-C — 1 process are obtained from Ginibre's theorem 14. 1. II 
by letting n — > oo. So from the point of view of determinantal processes, our 
problems can be stated as finding a probabilistic meaning to Det-S 2 — a and Det- 
B-q. 



32 



Chapter 5 

Matrix analytic functions on the 
sphere 

In this section we prove Conjecture 14.3. ll stated at the end of Chapter |4j 
i.e., we show that the processes Det-S 2 -a arise as the singular points of the matrix 
GAF zA — B or equivalently, zeros of the polygaf det(zA — B). Recall that for S7 = § 2 , 
a is a positive integer (the number of points in the process) . We shall denote it by 
n in this section. 

Theorem 5.0.4. Let A, B be independent nxn random matrices with Ltd. standard 
complex Gaussian entries. Then the set X of zeros of det(zA — B) = 0, i.e., the 
eigenvalues of A~ 1 B, has the distribution Det-§ 2 -n. 

We need the following lemma. 

Lemma 5.0.5. Let X be a point process on C with n points almost surely. Assume 
that the n-point correlation function (equivalently the density) of X has the form 

p(zi, ...,z n ) = \ A(zi, ...,z n ) | 2 l/(|zi| 2 , . . . , \z n \ 2 ). 

Here A(zi, ... ,z n ) denotes the Vandermonde factor f\ (zj — zA. 

i<j 

Suppose also that X has a distribution invariant under automorphisms of 
the sphere S 2 , i.e., under the transformations <p a ,p(z) = a -^ +l3 _ , for any a, (3 satisfying 
| a | 2 + |/3| 2 = 1. Then 

n 

F(M 2 , . . . , W 2 ) = Const. [J 7TX7-15VHT - (5.0. 1) 
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Proof of Lemma \5.0. 51 The claim is that the probability density of the n points of 
X (in exchangeable random order) is 

n 1 

q(z u ...,z n ):= Const. | A( Zl , . . . , z n ) \ 2 JJ 2 



fc=1 v i bi y 

First let us check that the density q is invariant under the isometries of § 2 . For 
this let a,B(z) = with a, (3 satisfying \a\ 2 + \/3\ 2 = 1. Then, 

= : * . 2 - (5-0.2) 

{-pz + a)^ 

1 4- I z\ 2 

l + \<p{z)\ 2 = - — _ 1 ' . (5.0.3) 

| — /9z + a | 

z — w 

cp(z) - <p(w) = - (5.0.4) 
{—pz + a)(—pw + a) 

From I5.0.2I . 15T0~31 and 15.0.41 . it follows that 

n 

q {<p{zx), p(z n )) Y[ W{z k )\ 2 =q(zi,...,z n ), (5.0.5) 

k=l 

which shows the invariance of q. 

Invariance of X means that Va, (3 with |a| 2 + |/3| 2 = 1, and for every z\, . . . , z n , 

we have 

n 

p((p(zx), . . . ,cp(z n )) \(f'(z k )\ 2 =p{z\, ...,z n ). (5.0.6) 

k=l 

Set W(zi, ...,z n ) = p h'-' z "| . Then, from 15.0.61 and 15.0.51 . we get 

• W (zi, . . . , z n ) is a function of \zk\ 2 , 1 < k < n, only. 

• W ((f(zi), <p(z n )) = W(z 1 ,...,z n ) for every z 1} ...,z n . 

We claim that these two statements imply that W is a constant. To see this fix 

Zk = ri e e' ldk , 1 < k < n, such that T\ < r k for k > 2. Let a = J— ? , 3 = Then 

\/l+rf V 1+r i 

\a\ 2 + |/3 1 2 = 1 and so ip a ^ is an isometry of S 2 . From the above stated properties of 
W , we deduce, 

W(z 1 ,...,z n ) = W((p(zi),...,<p(z n )) 

z 2 - z\ z n - Z\ 



WO, 



w[o, 



1 + Z 2 Zi 1 + Z n Z\ 

^e* 6 * 2 — z\ 
1 + r 2 e i ^z 1 



r n e L " n - zi 
1 + r n e l6n Zi 
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Take z\ = 1 and 1 < < 1 + e. Then as 2 < k < n vary independently over [0, 2ir], 
the quantities I r fc e ' efc - 21 i var y over intervals ^h-r, ^-'-r ■ By our choice of rus, 
this means that 



W(0, t 2 ,.-.,t n ) = Constant 



1 

e, - 

e 



e is arbitrary, hence W(0, t 2 , ■ ■ ■ ,t n ) is constant. Therefore W(z\, 
This shows that p(zi, . . . , z n ) = Const.g(^i, . . . , z n ). 



, z n ) is constant. 

□ 



Proof of Theorem \5. 0. 41 Firstly we observe that the distribution of X is invariant 
under conformal automorphisms of S 2 . This is a direct consequence of Proposi- 
tion I^^^D with Q ="det" and f(z) = az — b. Still we give another simple direct 
proof. Let a, (3 be such that \a\ 2 + |/?| 2 = 1. Set 



C = aA + f3B 



and 



D = -PA + aB. 



Then C and D are i.i.d. matrices with i.i.d. standard complex Gaussian entries. 
Therefore the eigenvalue set of C~ 1 D is also distributed as X. 
Now X is the set of solutions to the equation 

det (zA -B) = 0, 

say X = {z\, . . . , z n }. By our observation X also has the same distribution as the 
set of solutions to 

= det (zC -D) = det ((za + 0)A - (-]3z + a)B) , 

which is precisely {f a ^{zk)}- This proves the invariance. 

Now by Lemma I5.0.5L it suffices to show that the density of points in X 
is of the form given in 15.0. II . 

We use the following well known matrix decomposition. 
Schur decomposition: Any diagonalizable matrix M £ GL(n, C) can be written as 



M = U(Z + T)U*, 



(5.0.7) 



where U is unitary, T is strictly upper triangular and Z is diagonal. Moreover the 
decomposition is almost unique, in the following sense: 
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M = V(W + S)V* in addition to 15.0.71 . with V,S,W being respectively 
unitary, strictly upper triangular, and diagonal, if and only if the element of W 
are a permutation of the elements of Z, and if this permutation is identity, then 
V = UQ and QSQ* = T for some that is both diagonal and unitary, i.e., 6 = 

Diag(e iei , . . . ,e i6n ) . 



Corresponding to this matrix decomposition 15.0.71 , Ginibre 111) proved 
the measure decomposition 

Ginibre's measure decomposition: If M is decomposed as in 15.0.71 . with the 
elements of Z in a uniformly randomly chosen order, then 

Y[dm(M i:j ) = il[ \zi- Zj\ 2 Y[dm(z k ) J I IpmC^) I dv (U) (5.0.8) 

i,j \i<j k J \i<j J 

where v is a finite measure on the unitary group U(n) such that dv(UQ) = dv(U) 
for every diagonal unitary 6. 

Conditional on A, the matrix M := A^B has the density 

g -tr( M* A* AM) M|2n 

w.r.t. the Lebesgue measure on GL(n, C) C C . From the measure decomposition 
15.0.81 we get the density of Z, T, U, A to be 

J] \zi ~ Z,\ 2 H dm(z k ) j e -tr(A'MI+MM*)) | A | 2n -Q dm{Tij) "Q ^..^ 

v *<i k=l ) i<j i,j 

(We have omitted constants entirely) Thus the density of Z is obtained by inte- 
grating over T,U,A. Now write Z = QR where and R are diagonal matrices with 
the polar and radial parts of z^, respectively. Then 

MM* = UO(R + e*T)(R + e*T)*e*u*. 

As stated earlier, du(U@) = du(U). The elements of Q*T are the same as elements 
of T, but multiplied by complex numbers of absolute value 1. Hence, Q*T has 
the same "distribution" as T. Thus replacing U by Q*U and T by Q*T we see that 
the density of Z is of the form Y\i<j \ z i ~ z j\ 2 V(R)- This is the form of the density 
required to apply Lemma 15 .0.51 Thus we conclude that the eigenvalue density is 

n 

const, n i* - n ( i+ N 2)n+i - ( 5 -°- 9 ) 



i<3 
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To compute the constant, note that 

J n/n—l\ z k 

is an orthonormal set. Projection on the Hilbert space generated by these func- 
tions gives a determinantal process whose kernel is as given in the definition of 
Del - S 2 - n. Writing out the density shows that this is the same as the eigenvalue 
density that we have determined. Hence the constants must match. □ 
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Chapter 6 

Matrix analytic functions on the 
disk 

Coni ecture 14 . 3 . 2l at the end of Chapter|4]asserted that the singular points 
of the matrix GAF Aq + zA\ + z 2 A2 + ... or equivalently, the zeros of the polygaf 
det(^4o + zA\ + z 2 A2 + . . .) are distributed as Det — B — n. In other words they form 
a determinantal point process with kernel 

n 1 
K n {z,w) = — — T 

7T (1 — ZW) n+L 

w.r.t the reference measure dfi(z) = (1 — \z\ 2 )^ n ~ 1 ^ dm{z) on D. Here At are i.i.d. 
n x n matrices with i.i.d. standard complex Gaussian entries. 

As already emphasized, Proposition 12. 1. II applies to show that the zeros 
of the polygaf Aq + zA\ + z 2 A2 + . . . are stationary on the unit disk. In this chapter 
we shall prove that the first and second correlation functions agree with those of 
the determinantal process with kernel K n . 

First intensity: In the notation of Chapter^ we have Q = determinant, a homo- 
geneous polynomial in n 2 variables < n and fij(z) i.i.d. copies of the power 

oo 

series Yl a nZ n - 

n=0 

1 'roin 12.3.21 . the intensity of zeros is | L Alog(l — 1«| 2 ) -1 . By an elementary 
computation this comes out to be ^j^pp ■ Since this is the same as K n (z,z), it 
follows that the polygaf under consideration has the same intensity of zeros as 
the determinantal process with kernel K (we omit n in the subscript often). 
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2-point correlations: We shall prove that 

p 2 (z,w) - P i(z) Pl (w) = -\K(z,w)\\ (6.0.1) 

which shows that the 2-point correlations for the zeros of the polygaf agree with 
those of the determinantal process with kernel K. (Henceforth correlations are 
expressed w.r.t. the Lebesgue measure). 

We prove 16.0. li by using Theorem \5.0. 41 and a change of variables! 

Let gij(z) = aijZ + b^, where 0^,6^- are all i.i.d. CiV(0, 1). Then K(z,w) = 
1 + zw. Consider the polygaf G(z) = det(gjj(z)). This is precisely the polygaf 
considered in Chapter [BJ Thus we know from Theorem 15.0.41 (which we proved 
by certain matrix decompositions, not at all by using the formulas for correlation 
functions of polygafs) that the zeros of G are determinantal with kernel 

n (1 + zw) n ^ 1 

Mi + izi^a + M 2 )^' 

But the the general formula 12.3.31 for two-point correlations of zeros of 
polygafs applies to G also and hence it must be the case that 



(2?r) 2 'dzd^dwdzw ^ p (1 + zz)p(1 + ww)p ~ tt 2 (1 + zz) n+1 (1 + ww) n+1 ' 



(6.0.2) 

where C p depend only on Q, not on the GAFs that we feed in. 
From 16.0.21 we get 

(2tt) 2 9xi dy l dx 2 dy 2 ^ pl (1 + xm)P(l + x 2 m)P " vr 2 (1 + x 1 x 2 ^) n + 1 (1 + ym) n+1 ' 

(6.0.3) 

This is because, both sides of 16.0.31 are analytic in x\, x 2 and anti-analytic in yi,y 2 
and moreover, 16.0.21 says that the two are equal on the diagonal {x\ = yi,x 2 = y 2 }. 
Thus, by a standard (and elementary) fact that can be found in any intr oductory 



29), the two 



book on several variable complex analysis, see for example Rudin 
sides must be equal for all xi,x 2 ,yi,y 2 . 

Now make the substitution x\ = iz, x 2 = iw, y\ = iw, y 2 = iz to get 

16 d d d d ^ - 2 {l- zz)p{\-ww)p _ n 2 {I - zz) n ~ l {I - ww) n ~ l 



(27r) 2 dz&zdw&zw p^ 1 Pl (1 - zw)p(1 - wz)p n 2 (1 - zw) n+1 (l - wz) n+1 ' 



(6.0.4) 
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But again using 12.3.31 the left hand side is precisely what we get for p2(z,w) — 
pi(z)pi(w), when Q is the determinant of n 2 variables and fjj are i.i.d. copies of 

oo 

= Yl a nZ n , with a n being i.i.d. standard complex Gaussians. And the right 

n=0 

side of 16.0.41 is K(z, w)K(w, z). 

We already know that pi(z) = K(z, z). Therefore, the two point correlation 

p 2 (z,w) is 

K(z,z) K(z,w) 

det 

K(w,z) K(w,w) 
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Chapter 7 



Asymptotic Normality 



7.1 Background: Results for Gaussian analytic functions 

Sodin and Tsirelson 135) proved asymptotic normality for smooth (C^) 
statistics applied to the zeros of the three canonical models of Gaussian analytic 
functions in 11.2.91 . 11.2. 101 and 11.2. 1 II . as the density parameter L — ► oo. More 
precisely, they showed that for any real valued ip e C%(Q,), if 



then, 



and also that 



VVar(Z L (^)) 



Var(Z L (<^)) = £|| AVII^.j + o(2T 1 ), 

for a constant k that is described explicitly and the same for all the three geome- 
tries. (Note that as L — > oo, the variance goes to zero!) Here m* is the invariant 
measure on and A* is the invariant Laplacian. In other words 



dm*(z) 



dm(z) if fi = C 

t^L if^ = § 2 (7.1.1) 
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and 



A* 



A if O = C. 

(l + |z| 2 ) 2 A ifft = S 2 . 
(1 - H 2 ) 2 A if n = B. 



(7.1.2) 



7.2 Our results: for polygafs 

In this article we modify the method of Sodin and Tsirelson to obtain 
central limit theorems for smooth statistics of zeros of polygafs. One point of this 
exercise is to demonstrate that the random analytic functions can be studied at 
the level of generality introduced in Chap ter|IZ] (Recall that polygafs include matrix 
GAFs as very special cases) . 

Theorem 7.2.1. LetF L (z) = Q{i[ L \z), . . .,^ L \z)) where, 

• Q is a fixed non-random homogeneous polynomial in p complex variables, 

• f/ L) are Ltd. GAFs in 1 1.2.9V or 11.2.10V or 11.2. lit . 



ip : £1 — ► E is aCl function on Q. 



Set 



E Viz)- 
eP£ 1 {o} 



Then 



Moreover 



Z L (ip)-E[ZM] 
^Var(ZM) 

Var(Z L (cp)) = ^1||A>| 12 



N(Q, 1) 



+ o{L- 1 ) 



(7.2.1) 



£ ll— r llL2(m*) 

for a constant k(Q) that is described explicitly and is the same for all the three 
geometries. 

Let Zl((p) = Zl(^p) — J£[Zl{<p)]. Let x denote a standard (real) normal 
random variable. The idea of the proof is to show that 



E 



E 



Z L &fY 2 E[x s ]+E\z L (<p) 2 Y o(l) 



as L 



oo, 



(7.2.2) 
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for s = 1,2,3,.... Then the moments of Zh( }f} =JS^M^ converge to those of % an d 
convergence in distribution follows. To show 17.2.21 . we need to compute the 
moments of Zi((p). 

Recall the formula 12.2.21 



J (p(z)driF(z) = J A(p(z)-^- log \F(z)\dm('A 
n n 



(7.2.3) 



From this we can also write 



Z L (<p) = I A V (z)—log\F(z)\dm(z) 



where F(z) = ^^^ d/2 F(z) as defined in Section T2.31 Then one can write the mo- 
ments of Zl((p) as 



E^ L (^) s j = (2tt)- 
for s = 1,2,.... 



(7.2.4) 



i=i 



7.2.1 Central moments of 



From the homogeneity of Q, we can write Fl(z) = Q(g[ L \z), . . . , g^'^z)), 



where gi(z) 



J 1 ^ — . These are no longer analytic functions, but they are 
independent complex Gaussian processes on ft with constant variance 1 . Thus for 
any fixed z, we have that g\ L \z), 1 < i < p are i.i.d. standard complex Gaussians, 
and from 12.2.41 it follows that 



log|F L (*) 



E 



0, 



V m!n! 



rnk <z)^{z): 



k=l 



with coefficients C m n that are the same as in 12.2.41 and depend only on Q but 
not on z or even the GAF f . We get 



E 



n io gi^(^o 



!'».":, - V =1 V 3 3 



n E 

k=l 



3=1 



Here mo 



, m P j) and likewise for rij. 
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Each of the expectations on the right hand side of the above equation can 
be "evaluated" by the Feynman diagram formula 11.2.51 . We denote by Kl(z,w), 



the quantity E [g(z)g(io) 
of a Feynman diagram 
the Gaussian vector (g(^i), 



K(z,w) 



Also we write v(-/; z\, 



for the value 



of a Feynman diagram 7, with edge weights given by the covariance matrix of the 

,g(z s )). Then, we get 



E 



Y[\og\F L (zj) 

3=1 



e n 



c n 



{m j ,n i } 1 



<j<a 



m,!n,! , 



e n 



c n 



1 V m r n r 



Z\, . . . , z s j 



e n 



71v,7p \3 



v(U p k=l -f k ;z 1 ,...,z s ). 



Here in the last sum, 7^, 1 < k < p vary over all possible Feynman diagrams on 
labels {1,1, ... s,s} and nx, , rij G are such that the number of vertices labeled j 
in 7^ is mfcj and the number of vertices labeled j in j k is n;y. 

Put this together with 17.2.41 to deduce that E Z L (ip) s is equal to 




=T^T ) / f II A *M z i) ) w ( u Li7fc! z 1 ,...,z s )f[ dm(2 



(2k)- 

' — ' \ ""7 x/m,-!n,-! / J \ "T"7 " / """7 

(7.2.5) 

Here again, the sum is over all legal diagrams j k - 

Our goal is to prove 17.2.21 . To that end, we now consider the second 
moment, which is obtained by setting s = 2. We get 



E 



l[log\F L (z,) 

3=1 



e n 



c n 



molrii! 



v(v p k=1 ik;z 1 ,z 2 ), 



(7.2.6) 



7iv,7p \j=l \ ~~ 3' 3 ' 

where, each j k is now a Feynman diagram on vertices labeled {1, 1, 2, 2}. 

Now suppose s is even. Write 17.2.61 with z\,Z2 replaced by £2fc-i>22fc f° r 
k < § and multiply them together. On the right hand side, we get (the product 
of the values of Feynman diagrams is the value of the union of the Feynman 
diagrams) 

e (n 0m,n 



7l>— >7p 



w ( u fc=i7fc;^i, • • -,z s ), 
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where the sum is over all Feynman diagrams 7*. on vertices labeled by {1, 1, ... s, s}, 
such that, when all the vertices labeled by j,j are identified for each j, then U? =1 7fc 

has connected components {1, 2}, {3, 4}, . . . , {s — 1, s}. 

- s 
To get E [^(v?) 2 ] 2 , integrate against v^p- fl Azj<p(zj) w.r.t Lebesgue mea- 

s 

sure over s . This yields that E [2x(<£>) 2 ] 2 is equal to 

E (n 7== ) / \i[A Zj v(z j ))v(ul =l7k] z 1 ,...,z s )f[dm(z j ). 

(7.2.7) 

Here again the sum is over all diagrams that have connected components {1,2}, 
{3, 4}, . . . , {s — 1, s} (upon merging vertices labeled 

Instead of pairing {1, 2, ... , s} as {1, 2}, . . . , {s — 1, s}, we could use any other 
matching. Write the expression analogous to 17.2.71 for each matching, and add 
them all up. Recall that the number of matchings of {1, 2, ... , s} is E [x s ]. Thus we 
deduce that E [Z L (p) 2 ] % E [x s ] is equal to 



, ) / ( II A ^( z i) ) u ( u Li7fc; z 1 ,...,z s )f[ dm(zj] 



<*o- y. in;,;;; 

(7.2.8) 

where the sum is over all diagrams 7^, 1 < k < p such that Uj^ "splits", i.e., when 
vertices j,j are merged, we get s/2 component each of size 2. 

Compare 17.2.81 with 17.2.51 (the expressions are incomplete without the 
commentaries that follows after the equations!). The terms on the right hand side 
of 17.2.51 that are absent in 17.2.81 are precisely those, for which U^ =1 7fc does not 
split into s/2 components. The proof of 17.2.21 will be complete once we show that 

s 

these terms together contribute a negligible amount compared to E [Zl(</?) 2 ] 2 . 

7.2.2 Estimating the second moment 

This is the case s = 2, which we already dealt with in detail, in Chapter^ 
Particularly, from 12.3.31 we can write 



E 



i 1 00 r 

Z L {V?\ = 7^2 El^l 2 / ^)^{w)\k L {z,w)\ 2p dm(z)dm(w), (7.2.9) 
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where \C P \ 2 = V | C m ,n 1 2 - It is convenient to express everything in terms of 

m.=n.=p 

[he invariant quantities of f2. From 17. 1. II and 17. 1.21 we rewrite 17.2.91 as 



E 



\ ^ \C P \ 



A*^(z)A*ip(w)\K L (z,w)\ 2p dm*(z)dm*(w). 



Now split A>(z)AVH as (A*<p(z)) 2 + (A*<p(w)) 2 - (A*<p(z) - A*<p(w)) 2 . We get a 
sum of three integrals. The first two integrals are equal by symmetry. We shall 
argue that the last integral is negligible. Firstly note that 



\K L (z,w)\' 



if n = c 



|i+gjglff 

(T+lipHT+Rp) 



-v , i-u ,., t~t if ^ — S 2 



(7.2.10) 



(i-|*| 2 ) L (i-M 2 ) L 

ll-zu?! 2 ^ 



if n 



Thus \Ki(z, w)\ 2 is 1 when z = w and decays rapidly as (z, u>) moves away from the 
diagonal. Since (A*<p(z) — A*(p(w)) 2 vanishes on the diagonal, it is easy to calculate 
that 

' X (7.2.11) 



(A*(p(z) - A*ip(w)) 2 \K L (z,w)\ 2p dm*(z)dm*(w) = o — 

Lp 



The first two integrals give us 

(A*ip(z)) 2 \K L (z,w)\ 2p dm*(z)dm*(w). 



(7.2.12) 



Q 2 



Since is also invariant under isometries of f2, fixing z and integrating w.r.t w 
we get 

' \k L (0,w)\ 2p dm*(w) I I / (A*(p(z)) 2 dm*(z) I . 



Now from 17.2. 101 it can be checked by direct computation that 

if n = c. 



\K L (0,w] 



\2p 



M 



Lp 



(7.2.13) 



Lp-l 



if 



From 17.2.111 . 17.2.121 and 17.2.131 . we get 



E 



C n 



*,„l|2 



p=0 



lL 2 (m*) 



O — 



L 



(7.2.14) 
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This shows 17.2.11 with 



p=0 y 



7.2.3 Estimating the values of non-split diagrams 

We want to show that the contribution of non-split diagrams to ~E[Zl((p) s ] 
is negligible. Note that this includes all the diagrams in the case when s is odd. 
Consider any p-tuple of diagrams (71, ... , j p ) on labels {1,1, ... s,s} such that Likjk 
does not split into pairs when j,j are merged. We bound 

=l7fcj z li ■ ■ ■ , z s) JT dm(zj), 
I V =1 / i =1 

in absolute value by the obvious 

W^WL j \v(U p k=llk ;z 1 ,...,z s )\f[dm(z J ), 

i 3=1 

where A is the support of p. Then split the integral of ^(U^L^; zi,...,z s )\ as a 
product of integrals over the connected components of Uklk- 

Without loss of generality, let {1,2, ... ,r} be a connected component of 
Ufc7fc (when j, j are merged) . Since the weights of edges are bounded by 1 , deleting 
some edges will only increase the integral that we want to bound. We delete 
enough edges to get a spanning tree T on {1, 2, . . . , r}. So we are left with an 
integral of the form 

/ I"! \Kl(zi, Zj)\dm(zi) ...dm(z r ). 
At 

Integrate inwards starting with the leaves. From 17.2. 131 we get that 

/ II \^L(zi, Zj)\dm{z\) . . . dm(z r ) < CLT r . 

Multiplying the contribution from each component, we get 

/ "(uLi-r*; *i> •••,*«) n dm ^ < CL ~ s+no - of com P° nents of u ^. 

£ 3=1 

Non-split diagrams are precisely those that have less than | components, whence 

CP 1 

the right hand side is 0(L~ 



s-l 
2 
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Now we want to bound the total contribution of all unsplit diagrams. This 
can be done in the following manner. 

Approximate log |F(z)| by polynomials (by truncating the Wick expansion). 
For integration of A92 against a polynomial, our cruder bound on a single Feynman 
diagram suffices, since there are only finitely many terms and each of them goes 
to zero. So we get asymptotic normality for Aip integrated against polynomials. 
From this, one can deduce asymptotic normality for Zl(<p). We skip the details 
(see 35)). 



To summarize, we have argued that 



E 



E 



ZM 2 2 K[ X s ]=o(L 



-«/2\ 



and proved that E Z L (<p) 2 2 is of order L -s / 2 . Thus 17.2.21 follows. 



48 



Chapter 8 



Overcrowding Problems 



8.1 Statements of the problems 

In this chapter we go back to the canonical models of Gaussian analytic 
functions defined in Chapter^ Namely, consider the following Gaussian analytic 
functions (GAFs): 

• Planar GAF : The function defined in 11.2.91 . 



where as before a n are i.i.d. standard complex Gaussians. Almost surely, fx, 
is an analytic function in the unit disk (and no more) . 

We denote the zero set by Z. Let n(r) denote the number of points of Z in the disk 
of radius r around (The GAF will be clear from the context) . By the invariance 
of the zero sets, the results carry over to disks centered elsewhere. 

Yuval Peres asked the following question and conjectured that the prob- 
ability decays as e cm2l °g( m ) i n the planar case (personal communication). 




where a n are i.i.d. standard complex Gaussian random variables. 



• Hyperbolic GAFs : For each L > the function defined in 11.2. 1 II 




,71 
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Figure 8.1: Samples of the zero set of g. Left: The zero process sampled under 
certain sufficient conditions (see the conditions in the lower bound of the proof of 
Theorem 18. 1. It Take a = 0.5, r = 2,m = 16) on the coefficients forcing 16 zeros in 
the disk of radius 2. Right: The unconditioned zero process. 

Question: Fix r > 0, (r < 1 in the Hyperbolic case). Estimate P [n(r) > m] as 
m — > oo. 

One motivation for such a question is in Figure 15711 There one can see the 
distribution of the zero set under certain conditions on the coefficients that force 
large number of zeros in the disk of radius 2 (this is not the zero set conditioned 
to have overcrowding - that seems harder to simulate). The picture suggests that 
the distribution of the conditioned process may be worth studying on its own. A 
large deviation estimate of the kind we derive will presumably be a necessary step 
in such investigations. 

The answer is different in the two settings. We prove- 

Theorem 8.1.1. Consider the planar GAF g. For any e > 0, 3 a constant C2 (de- 
pending on e, r) such that for every m > 1, 

e -im 2 log(m)+0(m 2 ) < p[ n ( r ) > m ] < e)m 2 \og(m) ^ 

In particular, P[n(r) > m] = e -!™ 2 i°gM(i+°(i))_ 

Theorem 8.1.2. Fix L > and consider the GAFf^. For any fixed r < 1, there are 



50 



constants f3, C\,Ci (depending on L and r) such that for every m > 1, 

2 

d{r)e~W)\ < p[ n (r) > m] < C 2 (r) e - /3(r)m2 . 
We prove Theorem 18. 1. II in Section 18.21 Theorem 18. 1.21 in Section 18.31 

8.2 Overcrowding - The planar case 

In this section we prove Theorem 18. 1. II Before that we explain why one 
expects the constant | in the exponent in Theorem 18. 1. II by analogy with the 
Ginibre ensemble. 

8.2.1 Ginibre ensemble 

The Ginibre ensemble is the determinantal point process (earlier we de- 
noted this by Det-C — 1) in the plane with kernel 

K(z,w) = - e -¥ z \ 2 ~¥ w \ 2 + zW . (8.2.1) 

7T 

This process is of interest because it is the limit in distribution, as n — > oo, of 
the point process of eigenvalues of an n x n matrix with i.i.d. standard complex 
Gaussian entries (Theorem 14. 1. II . 

The Ginibre ensemble has many similarities to the zero set of g. In partic- 
ular, the Ginibre ensemble is invariant in distribution under Euclidean motions, 
has constant intensity ^ in the plane and has the same negative correlations as 
iJ g at short distances. Therefore there are other similarities too, for instance, 
see 0)- There are also differences between the two point processes. For instance, 
the Ginibre ensemble has all correlations negative, whereas for the zero set of g, 
long-range two-point correlations are positive. However, in our problem, since 
we are considering a fixed disk and looking at the event of having an excess of 
zeros in it, it seems reasonable to expect the same behaviour for both these point 
processes, since it is the short range interaction that is relevant. In case of the 
Ginibre ensemble, the overcrowding problem is easy to solve. 

Theorem 8.2.1. Let nc{r) be the number of points of the Ginibre ensemble in the 
disk of radius r around (by translation invariance, the same is true for any disk 
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of radius r). Then for a fixed r > 0, 

P [n<j(r) > m] = e7 



im 2 log(m) (l+o(l)) 



Proof. By Kostlan 21), the set of absolute values of the points of the Ginibre 
ensemble has the same distribution as the set {R±, R2, . . .}, where R n are indepen- 
dent, and R 2 n has Gamma(n, 1) distribution for every n. Hence R\ = £1 + . . . + £ n , 
where are i.i.d. Exponential random variables with mean 1, and it follows that 



P [Rl < r 2 ] > J] P 



k=l 



& < 



n 



,2 \ n 

> 'k ■ 



as long as n > r 2 , because P [£1 < x] > | for x < 1. Therefore we get 



P [n G (r) >m] > JJ P [i? 2 < 

n=l 
m 

a n 



n=l 



2 \ n 



2ra i 



m(m+l) m 

-2 \ 2 



E nlog(n) 

e n=1 



(8.2.2) 
(8.2.3) 

(8.2.4) 



Here and elsewhere we shall encounter the term Yl wlog(n). We compute its 

n=l 

asymptotics now. 



n log(n) < x log(a;) < (n + 1) log(n + 1) 



for n < x < n + 1 



Integrate from 1 to m + 1 and note that 

a 

y xlog(x)cte = ^a 2 log(a) - ^- + i, 
1 

to get 

E .. 1. \9 1 / \ (in -\~ l) 2 1 v — ^ , / \ 
nlog(n) < -(m + 1) log(m + 1) 1 — < wlog(n). 



n=l 

Thus 18.2.41 gives 



(8.2.5) 



n=l 



P[n G (r)>m] > e -*(™+U a lo g(m+ i)+^-± + ^±ii log(r 2 /2) 



e 2 



-m 2 log (m)+0(m 2 ) 
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To prove the inequality in the other direction, note that 



P [n G (r) >m]<P 



71=1 

2 



+ P i R n< 



r 21 



(2\ 
™)npK< r 2 ] + E'- 



Y,HRl<r 2 )>m 

n=m 2 +l 

-nlog(n)(l+o(l)) 

n=l n>m 2 

In the second line, for the first summand we used the fact that R\ are stochasti- 
cally increasing and for the second term we used the well known fact P [i? 2 < r 2 ] = 
P [Pois(r 2 ) > n] and then the usual bound on the tail of a Poisson random vari- 
able, namely P [Poisson((9) > a] < e -^o g (a/9)+a-e _ 

Using the same idea to bound P [i? 2 < r 2 ] in the first summand, we obtain 

m 

n\og(n/r 2 )—r 2 +n , p —m? log(m 2 )(l+o(l)) 



(2\ m 
m 1 Yl e -nlog(n/r2)-r 2 +n + g 
/ r, — 1 



n=l 

2 \ m(m+l) (r2)) _ mr 2_ £ nl (n) 

< I |e 2 n=l + e -m 2 log(m i )(l+o(l)) 



im 2 log(m)(l+o(l)) 



e 2 



(using 18.2.51 again) . 



In the last line we used (™ ) < m 2m . This completes the proof. □ 

8.2.2 Proof of Theorem 18.1.11 

Our method of proof is largely based on that of Sodin and Tsirelson 36). 
(They estimate the "hole probability", P [n(r) = 0] as r — > oo.) 

Proof of Theorem \8. l.ll Lower Bound Suppose the rrfi 1 term dominates the sum 
of all the other terms on dD(0;r), i.e., suppose 



— | > | - | whenever \z\ = r. (8.2.6) 

\im\ ^* vn! 



Then, by Rouche's theorem g(z) and a ^—^ have the same number of zeros in 
D(0; r). Hence n(r) = m. Now we want to find a lower bound for the probability of 
the event in 18.2.61 . Note that the left side of 18.2.61 is identically equal to 
Now suppose the following happen- 

1. \a n \ < n Vn > m+ 1. 
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2. \a m \ > (a + l)m where a will be chosen shortly. 

3. |au|-T=j < j t= for every < n < m — 1. 

1 1 vi' vrn! 

Then the right hand side of 18.2.61 is bounded by 

771—1 „ OO I 



RHSofdUl < V|a n |^=+ V 

n=0 n— ™ 1 



n=m+l 

oo ~ 

nr 



n=0 

r 



< m 



+ C 



n=m+l 

mr m 



ml ymi 
r 



(C + l)m 



'ml 



< \a r 



if a = C. Thus if the above three events occur with a = C, then the m^ 1 term 
dominates the sum of all the other terms on dD(0;r). Also these events have 
probabilities as follows. 

00 2 2 

1. P[KI <nVn>m + l] > 1- £ e" n > 1 - C'e"" 1 . 

n=m+l 

2. P[|a m | > (C+ l)m] = e -(^+ 1 ) 2 ™ 2 . 

3. The third event has probability as follows. Recall again that P [£ < x] > § 
if x < 1 and £ has exponential with mean 1. We apply this below with x = 
^ r m ^VnJ. ^j This i s dearly less than 1 if n > r 2 . Therefore if m is sufficiently 
large it is easy to see that for all < n < m — 1 , the same is valid. Thus 



r m-n./ n \ 



k < 



Vn < m — 1 



m— 1 



> 



n=0 
m— 1 

n 

n=0 



On I < 



n p 

r 2m— 2n^[ 



„m-n. /„| 



2m! 



„m(m+l)g^m 2 log(m) q— m^— in 2 Iog(m)+O(m 2 ) 



g— |m 2 log(m)+0(m 2 ) 



Since these three events are independent, we get the lower bound in the theorem. 
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Upper Bound By Jensen's formula, for any R > r we have 

R 2tt 2tt 

n(01og(f) < J^du = J log|g(^)|^-|log|g(re^)|^. 

r 



Let i? = i? m = y^m. Sodin and Tsirelson 36) show that 



loft .!/(/) .1' ( i + e ) t 2 



< e" 



(8.2.7) 



(8.2.8) 



where M(t) = max{|g(z)| : \z\ < t}. 

Now suppose n(r) > m and log M(R m ) < (| + e) m for some e > 0. Then by 
18.2.71 we have 

2tt 



y log |g(re* e )|^ > mlog ( - (\ + r ] //, 



-m log(m) — m log(r) — ( - + e ] m 



1 



mlog(m) — 0(m). 



Thus P[n(r) > m] is bounded by 



logM{R m ) > [ -+e)m 



+ P 



2- 



y log|g(re ie )|^ > -mlog(m) - 0(m) 



_ e _e + P 



2tt 



j log\g(re i9 )\^ > imlog(m)(l + o(l)) 



bv l8~2~71 



From Lemma 18.2.21 we deduce that for any 5 > 0, there is a constant C2 such that 

2tt 



1 

log|g(re ?e )|— > -mlog(m)(l + o(l)) 

Z7T Z 



< C 2e -(2-5)(f log(m)) 2 / log(f log(m)) 



< c , 2e" ( 5-|) m2lo s( m )( 1 +°( 1 )). 



From this, the upper bound follows. 



□ 



Lemma 8.2.2. For any given 5 > 0, 3C2 such that for every m, 



2tt 



log|g(re^)|— > m 
zir 



d6 



< C2& lo s( m ) . 
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Proof. Let P be the Poisson kernel on D(0; r). Fix e > and let A e = sup{P(re ld , w) : 
\ w \ = e ,e G [0,2vr)} and B € = mi{P(re i0 ,w) : \w\ = e,6 e [0,2vr)}. Since log |g| is a 
subharmonic function, for any w with \w\ = e, we get 

log|g(re^)|P(re^,u;) — 



2- 



_ -I I lo S _ |gK e )|^ - S 6 1 log_ |g(re^)|^ 



2tt 



2tt 



2tt 



2- 



d0 



•I / \o^\ S (n"')\^-+B e ( j'log\g(re ie )\^- /"log+ |g(re\, _ 



o 

2tt 



\0 







< 5 e / l og |g(re iy )|— +A £ log+M(r). 



2tt 



This implies logM(e) < 5 e / log |g(re ie )|ff + A e log + M( 



2zr 



Therefore if j log |g(re l(9 )||^ < —m, then one of the following must happen. 

o 71 

Either {logM(e) < -B t m + ^/m} or log + M(r) > y 7 ™}- 

Using 18.2.81 . since M(r) < M ICraJ for any C, we see that 



log + M(r) > 



< e" 



for some constant C depending on e. Hence 



2tt 



log|g(re ) — < -m 

Z7T 



< e' e ° m + P [log M(e) < -B t m + ^m] 



< g -)- g e log(m) 



r(l+o(l)) 



where in the last line we have used Lemma I8.2.3L 

As e —> 0, B e — > 1 and hence the proof is complete. 



□ 



Now we prove the upper bound on the maximum modulus in a disk of 
radius r that was used in the last part of the proof of Lemma I8.2.2L For possible 
future use we prove a lower bound too. 
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Lemma 8.2.3. Fix r > 0. There are constants a, C±, C 2 such that 

C ie iog(™o < p[l g M(r) < -m] < C 2 e~ I ^T ( + 1 Jj . 



fk-\ 

' I |2 

Lin 



1/2 



■ » . ~ |q n | r " 



+ E 



Proof. Lower bound By Cauchy-Schwarz, M(r) < I J2 

\n=0 / n=fc 

We shall choose /c later. We will bound from below the probability that each of 
these summands is less than 

Let ifk denote the density of T(k, 1). 

'k-1 \ V2 



e - 2 /2 < 1 



vn=0 



'fe-1 



e -2m g r 2 



.n=0 



-2m„r 2 \ „-2m„r 2 



e -2mfc-A;log(A;)+0(fc) 



Also if \a n \ < n 2 Vn > k, then the second summand 

r k k 2 



Vn7 ~ \/M " 



Also the event {|a„| < n 2 Vn > A;} has probability at least 1 - ^ e~ n > 1 — Ce~ k . 

n=k+l 

Thus if we set k = lo y g ^ for a sufficiently large 7, then both the terms are 
less than e~f with probability at least e -27™ 2 /i°sM. 
Upper bound By Cauchy's theorem, 

Vn! fg(C) 



2m J C n+1 



where C r is the curve C r (t) = re lt , < t < 2ir. Therefore, 



M(r)Vn\ 



Thus we get 



P[M(r) < e" m ] < J] P 



n=0 



k < 



- m Vn! 



\a n \ 2 are i.i.d. exponential random variables with mean 1. Therefore, 



K < 



< 



-2m n j 
r 2n 
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Using this bound for n < k := , , we get 



log(m) 

P[M(r) < e~ m ] < Yl 



I e -2m n | 



ra=0 

< Ce -2m/:+( log(fc)+0(fc 2 ) 



-2/? + is minimized when /3 = 2 and we get, 

P[M(r) < e" m ] < e- 2 ^^^^ 1 )). (8.2.9) 



□ 



8.3 Overcrowding - The hyperbolic case 

8.3.1 The determinantal case 

We give a quick proof of Theorem 18. 1 .21 in the special case L = 1, as it is 
much easier and moreover we get matching upper and lower bounds. The proof 
is similar to the case of the Ginibre ensemble dealt with in Theorem 18.2.11 and is 
based on the fact that the set of absolute values of the zeros of fi is distributed 
the same as a certain set of independent random variables. The reason for this 
similarity between the two cases owes to the fact that both of them are determi- 
nantal. The zero set of fi is a determinantal process with the Bergman kernel for 
the unit disk, namely 

K B (z,w) = -- 

. . IT (I — ZW) Z 

as discovered by Peres and Virag 128). 

Proof of Theorem \8. 1 .21 fo r L = 1. By the result of Peres and Virag quoted in Theo- 
rem 14.1 .21 the set of absolute values of the zeros of fi has the same distribution 
as the set {Un 2n } where U n are i.i.d. uniform[0, 1] random variables. Therefore, 

m 

/2n 



P [n(r) > m] > [J P \ul 



< r 

n=l 



r 2n 



n 

n=l 
m(m+l) 
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To prove the inequality in the other direction, note that 



P [n(r) > m] < P 



< 



m 
m 



£ l{U l J 2n <r)>m 

m 

n p K /2n 



n=l 
2 



OC 



+ £ P t/V 2 " < 



< r 



n=m 2 +l 
In 



+ E r 

n>m 2 



> m(m+l) 

m 



+ 



l_ r 2 

_ r m(m+l) ^1 _|_ O ^ e ml °g( m )^ _ 

This completes the proof of the theorem for L = 1. 



□ 



8.3.2 All values of L 



Remark 8.3.1. Overall, the idea of proof is the same as in that of Theorem 18. 1 . Tl 
However we do not get matching upper and lower bounds in the present case, the 
reason being that in the hyperbolic analogue of Lemma 18.2.31 the leading term 
in the exponent of the upper bound does depend on r, unlike in the planar case. 
(An examination of the proof of Theorem 18. 1. II reveals that we get a matching 
upper bound only because replacing r by e does not affect the leading term in the 
exponent in the upper bound in Lemma 18.2.31 . However we still expect that the 
lower bound in Theorem 18. 1.21 is tight. (See remark after the proof). 



Proof of Theorem \8. 1 .21 Lower Bound As before we find a lower bound for the 
probability that the m term dominates the rest. Note that if \z\ = r, 



-V 

rn 



n1/2 m-l /_r\l/2 00 /-A 1/2 

I a m z m \ < KM ) r n + Y, Ml n ) rn ( 8 - 3 -^ 

' n=0 ^ ' n=m+l ^ ' 



Now suppose the following happen- 

1. \a n \ < \fn Vn > m + 1. 

2. \a m \ > (a + l)^Jm where a will be chosen shortly. 

3 - \ a n\{'n) 1/2rn < ^m(~m) 1/2rm for evei y < n < m - 1. 
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Then the right hand side of 18.3. 11 is bounded by 

m - 1 /_r\ 1 / 2 00 /-A 

RHSofjasU < eki L rn + E ki n 



n=0 
m— 1 



< 



n=0 



-L 



m \ m 



-L 

m 



1/2 



n=m+l 
1/2 oo 

r m + E v " 

n=m+l 



ii 



1/2 



1/2 



(C + l)xM 



r m + CVm" 



r ' 



for some C 



m 



< \a r 



-A 



1/2 



if a = C. Thus if the above three events occur with a = C, then the rrfl 1 term 
dominates the sum of all the other terms on dD(0;r). Also these events have 
probabilities as follows. 



1. P[M < -v/n Vn > m + 1] > 1- £ er n >\-C'e 

n=m+l 



2. P[|a m | > (a + l)vH =e 



-(a+l) 2 m 



3. The third event has probability as follows. Recall again that P [£ < x] > | if 
x < 1 and £ has exponential distribution with mean 1. We apply this below 



with x 



-L\l/2 



(if) 



1 | < Km> - 1/2 Vn < m - 1 



. This is clearly less than 1 . Thus 

V m } 



m 



m— 1 



> 



rip 

n=0 

71—1 

n 



L\ 1/2 ^m— n 



On < 



ml 



-Ml/2 



m-1 c-M 2m-2n 



n=0 



2 ™(n) 



> r 



m(m+l) m -m 



m—1 

n 

n=0 
m—1 

n 



(m + 1) . . . (m + L - 1) 
(n+l)...(n + L-l) 



m 



n=0 

-> „m(m+l)+0(mlog(m)) 



(n + L) L 



Since these three events are independent, we get the lower bound in the theorem. 
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Upper Bound The proof will proceed along the same lines as in Theorem 18. 1. 11 
We need the following analogue of Lemma 18.2.31 

Lemma 8.3.2. Fixr < 1. LetM(r) = sup^ eD(0;r) Then 

P[M(r) < e~ m ] < e ~WJJ {1+o{1)) . 
Proof. By Cauchy's theorem, for every n > 0, 



n 



2m J C n+1 

rT 



From this we get 



Since ( n L ) > r(L+i) 1 we orj tain 



(n) 



P[M(r) < m] < [JP[|a n | 2 < 



2 r(L + l)e" 



-2m 



< 



log(l/r) 

n 

n=0 



r(L + l)e~ 2m 



< e 



-I5i^ + (l^TM) 2l °g(r)+0(mlog(m)) 

- logTl 2 /r) +0( " ll0g( " l)) . 



□ 

Coming back to the proof of the upper bound in the theorem, fix R such 
that r < R < 1. Then by Jensen's formula, 



n(r)]og(!) < | ^d«= | log|/(^)|£-yiog|/(re^)|^. 

r RT rT 

Now consider the first summand in the right hand side of 18.2.71 . 



(8.3.2) 



jev de 



io g |f(fte w )|— > 



LRT 



< P [log M(R) > 



rn 
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Now suppose that \a n \ < \ n Vn > m + 1 where 1 < A < 1/R. This has probability at 
least de- x2m l 2 . Then, 

1/2 



M(R) < ^2 KU n J R 



n-0 
m 



1/2 



\n=0 



for some constants Cr and Cr>. 



Thus if M(i?) > ev 7 ™ then either £ K\ 2 > Ce 2 ^™ or else |a n | > A n for some 



n > m + 1. Thus 



This proves that 



n=0 



P [M > v 7 ^] < e" 



\og\i{Re ie )\— > 

Z7T 



LRT 



< e" 



Fix 5 > and close enough to 1 such that log(-R) > — <5. Then with probability 
> 1 - e- £cV ^, we obtain from 18.3.21 . 



,d0 



1 2tt 



log \f(re l0 )\— > m log - - S - y 7 ^- 



Now the calculations in the proof of Lemma 18.2.21 show that 



2tt 



log |f (re id )\P(re ie , w)— + A e log, M(r). 

27T 



Here < e < r is arbitrary and A e , B e are as defined in Lemma T8.2.2L By the same 
computations as in that Lemma, we obtain, we obtain the inequality 

2tt 



/ 



d9 

log |f(re* 9 )| — < -m(| logr| - 5) + \fm 



< e ^ +e -e _ 



Therefore, by 18.3.21 



P [n(r) > m] < e 



-Km 2 log 2 (r)(l+o(l)) 



where k = sup 



Br 



< e < r > . However it is clear that this cannot be made to 



match the lower bound by any choice of e. 



□ 
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Remark 8.3.3. If we could prove 

that would have given us a matching upper bound. Now, one way for the event 

J log \i(re ld )\^ < —x to occur is to have log M(r) < —x which, by Lemma 18.3.21 

o 

has probability at most e~ x / lo s(?). One way to proceed could be to show that 
if the integral is smaller than — x, so is logM(s) for s arbitrarily close to r (with 
high probability) . Alternately, if we could bound the coefficients directly by the 
bound on the integral (as in Lemma 18.3.21 , that would also give us the desired 
bound. For these reasons, and keeping in mind the case L = 1, where we do have 
a matching upper bound, we believe that the lower bound in Theorem 18. 1.21 is 
tight. 



2tt 



|log|f(re*)||U 
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Chapter 9 

Moderate and very large 
deviations for zeros of the planar 
GAF 



Inspired by the results obtained (using not entirely rigorous physical ar- 
guments) by Jancovici, Lebowitz and Manificat for Coulomb gases in the 
plane (eg., Ginibre ensemble), M.Sodin 1341) conjectured the following. 

Conjecture 9.0.4 (Sodin). Letn(r) he the number of zeroes of the planar GAFg in 
the disk D(0, r). Then, asr^oo 



log log 



P[\n(r)-r' 2 \>r a ] 



log r 



i 



< a < 1; 



2a — 1, 2 
3a - 2, 1 < a < 2; 
2a, 2 < a. 



(9.0.1) 



The idea here is that the deviation probabilities undergo a qualitative 
change in behaviour when the deviation under consideration becomes comparable 
to the perimeter (a = 1) or to the area (a = 2) of the domain. 

Sodin and Tsirelson 13611 had already settled the case a = 2 by showing 
that for any 5 > 0, 3ci(<5), 02(6) such that 

e -ci(5)r 4 < p[| n ( r ) - r 2 \ > 5r 2 } < e~ C2 ^ r \ 
Here we consider P[n(r) — r 2 > r a ] and prove that a "phase transition" in 
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the exponent occurs at a = 2. More precisely we prove that the conjecture holds 
for a > 2 and show the lower bound for 1 < a < 2. 

Theorem 9.0.5. Fix a > 2. Then 

P [ n (r)>r 2 + 1 r a ] = e -(f-i)7 2 ^ lo g r(i+o(i))_ 

Theorem 9.0.6. Fix 1 < a < 2. Then for any 7 > 0, 

P [ n{r ) > r 2 + 7r a] > e - 7 3 r 3«-2 (1+o(1)) _ 

We prove Theorem l9~05l in Section 19771 and Theorem l9~076l in Section l9~2l 
Taken together these show that the asymptotics of P \n(r) > r 2 + jr a ] does un- 
dergo a qualitative change at a = 2. 

Remark 9.0.7. Nazarov, Sodin and Volberg have recently proved all the remaining 
parts of the conjecture (personal communication). 



9. 1 Very large deviations for the planar GAF 

In this section we prove Theorem l9.0.51 

Remark 9.1.1. In the case a > 2, one side of the estimate as asked for in the 
conjecture (with log log of the probability) follows trivially from the results in Sodin 
and Tsirelson l36D . They prove that for any 5 > 0, there exists a constant c(5) such 
that 

P [\n(r)-r 2 \ > 5r 2 ] < e -< S)r \ 

When a > 2, clearly n((l — 5)r^) > n(r), whence from the above result it follows 
that 



P [n(r) >r 2 + r a ] < P [n((l - 5)r^*) > r c 

-c{6)r 2a 



< e" 



This gives 



log log ( p[| n ( r )^ r ^| >r . Q ] 

limsup — < 2a. 

r^oo log r 



(9.1.1) 
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The obviously loose inequality n((l — 8)r^) > n(r) that we used, suggests 
that 19. 1. II can be improved when a > 2 to Theorem 19.0.51 



(hnz | > | ^"fL | whenever \z\ = r. (9.1.2) 

\fm\ vn! 



Proof of Theorem \9.0.5[ Lower Bound Let m = r 2 + jr a . Suppose the irfi 1 term 
dominates the sum of all the other terms on dD(0;r), i.e., suppose 

(I ,,, l" , ^-r fflri, 

Now we want to find a lower bound for the probability of the event in 
19.1.21 . Note that the left side of 19.1.21 is identically equal to 
Now suppose the following happen- 

1. | On, | < n Vn > m + 1. 

2. \a m \ > m. 

3. \a n \ < lllj^ f or every < n < m — 1. 

Vn! m vm! 

Then the right hand side of 19. 1.21 is bounded by 

m— 1 n oo I I n 

RHSofeXH < V|a n |^=+ V 

vn! — , , vn! 

n=0 n=m+l v 

E7r r v — v nr 
i ^+ 7 -r^ 
„ m Vm! ! Vn! 

n=U ra=m+l v 

^ mr m / jr a 
~ Tn^I V~ ° 

^ Qm r 
m! 

Thus if the above three events occur, then the rrfi 1 term dominates the 
sum of all the other terms on dD(0;r). Also these events have probabilities as 
follows. 

00 2 2 

1. P[\a n \ <nVn>m + l] > 1- £ e~ n > 1 - C'e"" 1 = 1 - o(l). 

n=m+l 

2. P[|a m | > m] = e~ m2 = e -7 2 ^(i+o(i)). 

3. The third event has probability as follows. Recall again that P [£ < x] > § 
if x < 1 and £ is exponential with mean 1. We apply this below with x = 
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7 2 f 2a r 2m-2ri n | 



This is clearly less than 1 if n > r 2 . Therefore if m is sufficiently 



large it is easy to see that for all < n < m — 1, the same is valid. Thus 



\a n \ < ^= — Vre < m — 1 

m 



m—l 



n=0 



On. < 



m ~ 1 J2 r 2a r 2m-2n | 
II 



n=0 



2m! 



r 2a(m+l)+m(m+l)2-ni m -2m e £ 



m 2 log(r)— hm 2 log(m)+0(m 2 ) 



e -(f -l)l 2 r 2a log(r)+0(r 2a ) 



Since these three events are independent, we get 

p [„( r ) > r 2 + r aj > e -(f-l) 7 V«logr+0(r 2 «)_ 



(9.1.3) 



Upper Bound We omit the proof of the upper bound, as it follows the same lines 
as that of Theorem 18. 1. II and we have already seen such arguments again in the 
proof of Theorem 18. 1.21 (In those two cases as well as the present case, we are 
looking at very large deviations, and that is the reason why the same tricks work). 

Moreover note that the lower bound along with 19. 1. II proves the state- 
ment in the conjecture. □ 



9.2 Moderate deviations for the planar GAF 

In this section we prove Theorem 19.0.61 

Proof of Theorem \9. 0. 61 Write m = r 2 + jr a . As usual, we bound P [n(r) > m] from 
below by the probability of the event that the m term dominates the rest of the 
series. 

Firstly, we need a couple of estimates. Consider ^ as a function of n. 
This increases monotonically up to n = r 2 and then decreases monotonically. 
m = r 2 + jr a is on the latter part. Write M = r 2 — jr a . 
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Firstly, observe that (r 2 — k)(r 2 + k) < (r ) , for 1 < k < ^r a , whence 

m— 1 

r 2m-2M > Y\ j. This implies that 

j=M+l 



VM! Vm\ 
Secondly, note that for any n = M — p, 



<^= (9.2.1) 



r 2n /n\ P t4 M -j 

n 



r 2M /M\ 11 r 2 

7 ' i=o 



p-1 

= \{{l-ir a - 2 -jr- 2 ) 

3=0 
p-i 

- E {yr a ~ 2 +jr- 2 ) 
< e J=0 



_ 7pr . Q -2_£(i^i) r .-2 



Now we set p = Cr 2 a with C so large that e 7C < \ . 

Then also note that if n < M — kp, it follows that 



Jin j | i 

7 <i, 0-2.2) 



r 2m ^/ m ] — ^k ' 

where we used 19.2. II to replace M by m. 

Thirdly, if n = m + p with p <r 2 — jr a , then, 

r 2n /n\ ^ r 2 

r 2m /ml -f ■ —^m-\-j 

p 

-3 E (7r Q - 2 +J>- 2 ) 

< e J=1 

= e -|(7pr"- 2 + 2^r- 2 )_ 

If p = 2Cr 2 ~ a , where C was as chosen before, then for n > m + kp, we get 

r2 " /n] < i. (9.2.3) 



r 2m / m \ — ^k 

From now on p = 2Cr 2 ~ a is fixed so that 19.2.21 and 19.2.31 are satisfied. 

Next we divide the coefficients other than the rrfi 1 one into groups: 
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• A k = {n : n G (M - kp,M - (k - l)p] for 1 < k < \f~}}. 

• Dk = {n : n G [m + (k — l)p, m + kp) for 1 < k < py] }• 

• B = {n : n G [M + l,m - 1]}. 

• C = {n : n G [2r 2 , oo)}. 

Remark 9.2.1. As defined, there is an overlap between DrMi and C. This is 



inconsequential, but for definiteness, let us truncate the former interval at r 2 
(just as ArMi is understood to be truncated at 0). 



Now consider the following events. 

1. K| < £ forn G A k for k < \f] }. 

2. |a n | < g for n G L> fc for k < \f]}. 

4. |a„| < re — 2r 2 for re G C 

5. |a m | > 15. 

Suppose all these events occur. Then 
1. The event \a n \ < g for n G /c < |~— ] gives 




(9.2.4) 




by (19.2.21) (9.2.5) 




\M/p\ 



(9.2.6) 



(9.2.7) 
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2. The event \a n \ < jt for n G D k , k < \—] gives 



2r 2 „ \M/p] M+kp 

a n z n 



8 up { | E 7=ri ; N= r > = E E i-i^j t9- 2 -8) 

n=m+l v k=l n=m +(k-l)p+l V 

* £^17 by (9.2.9) 



7^ V + ^ 



< -^=fl + ■£:)• (9.2.10) 



3. The third event gives 



71 III 



by assumption. 
4. The event \a n \ < n — 2r 2 for n £ C: Since n > 2r 2 , 



< 



< 





n 

fc=m+l 


r 

71 


j,m 
\fm). 


n 

n 

fc=2r 2 + 


r 






n-2r 2 






1 



Therefore we get (using \a n \ < n — 2r Vre > 2r ) 



r n m / i \ n-2r 2 



neC v " v ' n>2r 2 

r m V2 



I(V2-i) 2 ' 



(9.2.13) 



Putting together the contributions from these four groups of terms, and 
using \a m \ > 15, we get (for large values of r) 



£ K\^= < \a 



r 



Now we compute the probabilities of the events enumerated above. 
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1. The event \a n \ < 4j for n £ A k for k < \^-~\. Now for a fixed k < 31og 2 (r), we 



deduce 



\a n \ < for n € A k 



> P 



> 



I «n | < for n € A k 
1 V 



2M 2 J " 



Therefore 
P 



a n \ < — for n £ A k for every k < 31og 2 (r) 



for some c. 

Next we deal with k > 31og 2 (r). 



> 



\ 3plog 2 (r) 



2M 2 



> e _ cr 2- Q(log(r)) 2^ 



I On I < — for n £ A k for every k > 31og 2 (r) 



2 An 

l a l > T7 
' ' M 



2 2k M- 



> 1- J] pP 

fc>31og 2 (r) 

= 1 - ^ Pe ~ 

fc>31og 2 (r) 

Now the summation in the the last line has rapidly decaying terms and starts 



with pe 



_ 2 61og 2 (r) Jlf-2 



which is smaller than pe r . Thus 



I On I <-gforn£4 for every k > 31og 2 (r) 



l-o(l). 



Thus the event in question has probability at least e CT " 2 Q ( 1 °s( r )) 2 ( 1 +°( 1 )). 



2. The event \a n \ < 4j for n e D fe for fc < [^] . Following exactly the same steps 



as above we can prove that 



I o n | <^forn£ D k 



> e -cr- 2 -«(log(r)) 2 (l+o(l))_ 



3. The event ^ |a„|-^= < 4^=. By Cauchy-Schwarz, 



rte_B 



E Kl 



< 



|2 \ / ^r 2 ™ 



\neB 



KneB 
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Y = £ Kl 2 has r(|5|, 1) distribution. Also £ ^ < e r , since the left hand 

2 

is part of the Taylor series of e r . Therefore the event in question has proba- 
bility, 



P [event in question] > P 



„2m 



Y < 16- 



„2m 



> 8 e -r 8 



?n! 



m! 



where 99 is the density of the r(|5|, 1) distribution. This last follows because (p 
is increasing on [0, \B\] and thus P[Y < x] > y> (§) f . for a; < \B\. Continuing, 



1 «r 2 ™ -, 2 / r 2 " 1 

Pfevent in question] > — ^ ^ [ 8 ^V e ~ r ) ^ 21 ^ 

^> ^<g27r a (mlog(r 2 )— r 2 — mlog(m)+m)+0(r a log(r))jg 2 JgJ 

where we used Stirling's approximation. 

The exponent needs simplification. Take the first and third terms in the 
exponent. We have —2jmr a log (™) . Recall that m = r 2 + jr a and that a < 2. 
Therefore by Taylor's expansion of log(l + jr a ~~ 2 ) we get 

(77 
— 
r 

Now consider 19.2. 151 . Expand the fourth term in the exponential as 2^r 2+a + 
2^ 2 r 2a . We get the following terms 

(a) r 2+a (—2j + 27) = 0, from the second and fourth terms (first piece of the 
fourth term) in the exponential in 19.2. 151 . 

(b) r 2a (— 27 2 + 27 2 ) = 0, from the sum of the first term in the expansion 
19.2. 161 and the second piece of the fourth term in the exponential in 
19.2.151 . 

(c) r 3a - 2 (^ 3 _ 27 s ) = _^,3 r 3a-2^ f rom y-^ expansion 19.2. 161 . 

(d) Other terms such as r a log(m), r a log(r), r a , r 4a_4 , r 5a_6 etc. All these are 
of lower order than r 3Q 2 when 1 < a < 2. 



_2 7 V a 
_2j 3 r 3a ~ 2 



_|__ / 3 r 3a-2 
+7 4 r 4a-4 



2 7 4 r 4a-4 + _ 



2^5^,50?— 6 _|_ 



(9.2.16) 
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Hence, 

P[event in quesion] > e -7 3 r- 3t *- 2 (i+o(i)). 

4. The event \a n \ < n — 2m for n £ C. This is just an event for a sequence of i.i.d. 
complex Gaussians. It has a fixed probability p (say). 

5. The event \a m \ > 15 also has a constant probability (not depending on r, that 
is). 

This completes the estimation of probabilities. Among these five events, the third 
one, namely V \a n \^Q < 4^= has the least probability (Recall that 1 < a < 2). 

Also these events are all independent, being dependent on disjoint sets of 
coefficients. Thus P [n(r) > r 2 + 7 r Q ] > e -7 3 r 3a - 2 (i+o(i))_ 

□ 
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